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LECTURE fiERTPK "BOUNDARY LAYER THEORY” 

Part H — Turbulent Flows* 

By H. Schlichting 

CHAPTER YTTT . GENERAL REMARKS ON TURBULENT FLOWS 
a. Turbulent Pipe Flow 


The flow laws of the actual flows at high Reynolds numbers differ 
considerably from those of the laminar flows treated in the preceding 
part. These actual flows show a special characteristic, denoted as 
"turbulence.” 

The character of a turbulent flow is most easily understood in 
the case of the pipe flow. Consider the flow through a straight pipe of 
circular cross section and with a smooth wall. For laminar flow each 
fluid particle moves with uniform velocity along a rectilinear path. 
Because of viscosity, the velocity of the particles near the wall is 
smaller than that of the particles at the center. In order to maintain 
the motion, a pressure decrease is required which, for l a mi nar flow, is 
proportional to the first power of the mean flow velocity (compare 
chapter I, Part I) . Actually, however, one obfeerveB that, for larger 
Reynolds numbers, the pressure drop increases almost with the square of 
the velocity and is very much larger than that given by the Hagen— 
Poiseuille law. One may conclude that the actual flow is very different 
from that of the Poiseuille flow. 

The following test, introduced by Reynolds, is very instructive: 

If one inserts into the flowing fluid a colored filament one can observe, 
for s-mal l Reynolds numbers, that the colored filament is maintained down- 
stream as a sharply defined thread. One may conclude that the fluid 
actually flows as required by the theory of laminar flow: a gliding 

along, side by side, of the adjoining layers without mutual mixing 
(laminar = layer flew) . For large Reynolds numbers, on the other hand, 
one can observe the colored filament, even at a small distance downstream 
from the inlet, distributed over the entire cross section, that is, mixed 


*"Vortragsre ihe 'Grenzschichttheorie. 1 Teil B: Turbulente Stromungen. ” 

Zentrale fur wissenschaftliches Berichtswesen der Luftfahrtforschung des 
Generalluftzeugmeisters (ZWB) Berlin-Adlershof , pp. 154—279* ’ The 
original language version of this report is divided into two main parts, 

Teil A and Teil B, which have been translated as separate NACA Technical 
Memorandums, Nos- 1217 and 1218, designated Part I and Part II, respectively. 
This report is a continuation of the lecture series presented in part I, the 
equations, figures, and tables being numbered in sequence from the first 
part of the report. For general information on the series, reference 
should be made to the preface and the introduction of Part I. 
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to a great extent with the rest of the fluid. Thus the flow character 
has changed completely for large Reynolds numbers: A pronounced 

transverse mixing of adjacent layers takes place. Irregular additional 
velocities in the longitudinal and transverse' directions- are superposed 
on the main velocity. This state of flow is called turbulent. As a 
consequence of the mixing the velocity is distributed over the cross 
section more uniformly for turbulent than for laminar flow (compare fig. 17 , 
part I) . For turbulent flow there exists a very steep velocity increase 
in the immediate neighborhood of the wall and almost constant velocity in 
the central regions. Consequently the wall shearing stress is considerably 
larger for turbulent than for laminar flow] the same applies to the drag. 
This follows also from the. fact that in turbulent flow a considerable 
part of the energy is used up in maintaining the turbulent mixing motion. 

The exact analysis of a turbulent flow shows that at a point fixed 
in space the velocity is subjected to strong irregular fluctuations with 
time (fig. 72) . If one measures the variation with time of a velocity 
component at a fixed point in space, one obtains, qualitatively, a 
variation as shown in figure 72. The flow is steady only on the average 
and may be interpreted as composed of a temporal mean value on which the 
irregular fluctuation velocities are superimposed. 


The first extensive experimental investigations were oarried out 
by Darcy (reference 60) in connection with the preliminary work for a 
large water-distributing system feu* the city of Paris. The first quanti- 
tative experiments concerning laminar pipe flow were made by Hagen 
(reference 95) • The first systematic tests regarding the transition from 
laminar to turbulent flow were made by Osborn Reynolds (reference 6l) . 

He determined by experiment the connection between flow volume and pressure 
drop for turbulent flow and investigated very thoroughly the transition 
of the laminar to the turbulent form of flow. He found, in tests of 
various velocities and in pipes of various diameters, that transition 

always occurred at the same value of the Reynolds number: This 

Reynolds number is called the critical Reynolds number. The measurements 
gave for the pipe flow: 


- (flrlt " 2300 


(13.1) 


For Re < Re cr ^ the flow is laminar, for Re > Re cr ^^., turbulent. Later 
on it was ascertained that the numerical value of Re cr it is, moreover, 
very dependent on the particular test conditions. If the entering flow 
was very free of disturbances, laminar flow could be maintained up to 
Re = 2^000. However, of main interest for the technical applications 
is the lowest critical Reynolds number existing for an arbitrary disturbance 
of the entering flow, due either to irregularities in the approaching flow 
or to "vortices for ming at the pipe inlet. Concerning thejlrag law of the 
pipe Reynolds found that the pressure drop is proportional to the 1.73 
power of the mean flow velocity: 
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1>. Turbulent Boundary Layers 


Recently It was determined that the flow along the surface of a body 
(boundary layer flow) also can be turbulent. We had found, for instance, 
for the flat plate in longitudinal flow that the drag for laminar flow iB 

proportional to \J U Q ^ (compare equation (9-18), Part I.) However , towing 
tests on plates for large Reynolds numbers carried out by Froude 

I.85 

(reference 96) resulted in a drag law according to which W ~ U 0 . More- 
over, the drag coefficients in these measurements remained considerably 
higher than the drag coefficient of the laminar plate flow according to 
equation (9. 19), Part I. Presumably this deviation is caused by the 
turbulence of the boundary layer. — 


A clear decision about the turbulent flow in the boundary layer was 
obtained by the classical experiments of Eiffel and Prandtl concerning 
the drag of spheres in 1914- (reference 62) . These tests gave the following 
results regarding the drag of spheres (compare fig. 73). The curve of 
drag against velocity shows a sudden drop at a definite velocity T cri ^-, 

If one plots 

(3? = frontal area) against the 
decrease to 2/5 of its original value 
) . Prandtl explained this phenomenon 


although it rises again with further increasing velocity 
the drag coefficient c^. = W/F — 

U Q d/v, Cv shows a 


Reynolds number U 0 d/V, c^ 
at a definite Reynolds number 


^ e crit 1 

in 191^. He was able to show that this drag decrease stems from the laminar 
boundary layer changing to turbulent ahead of the separation point. The 
resulting considerable rearward shift of the separation point causes a 
reduction of the vortex region (dead water) behind the sphere (fig. 74). 

This hypothesis could be confirmed by experiment: by putting a wire ring 

on the sphere (sphere diameter 28 centimeters, wire diameter 1 millimeter) 
one could attain the smaller drag at smaller V cr and Re cr if The wire 

ring is put on slightly ahead of the laminar separation point; it causes a 
vortex formation in the boundary layer, which is thus made turbulent ahead 
of the separation point and separates only farther toward the rear. By 
meanB of the wire ring the boundary layer is, so— to— say, "infected" with 
turbulence. Due' to the mixing motions which continually lead high velocity 
air masses from the outside to the wall, the turbulent boundary layer is 
able to overcome, without separation, a larger pressure increase than the 
laminar boundary layer. 


The turbulence of the friction layer is of great importance for all 
flows along solid walls with pressure Increase (diffuser, wing suction 
side). It is, however, also present in the flow along a flat plate where 
the pressure gradient is zero. There the flow in the boundary layer is 
laminar toward the front, experiencing transition to the turbulent state., 
further downstream. Whereas the laminar boundary layer thickness increases 
downstream with xl/2, the turbulent boundary layer thickness increases 
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approximately as x^; that is, for the turbulent boundary layer the 
increase of the boundary layer thickness is considerably larger, (fig. 75 ). 
The position of the transition point x^^^. is given by (fig. 75): 


( 


= 3 to 5 x'lo5 

'crit 


(13-2) 


Hi comparing the critical Reynolds numbers for the pipe and the plate 
one must select r and S, respectively, as reference lengths. The 
equation for the flat plate is, according to Blasius (reference 8 ) 
(compare equation ( 9 * 21 a)) 


or 


5 = 


M - 5 


= 5\ — 



IU 0 x 


Thus, with 


(¥> 


(13-3) 


= 25000 : 


crit 

m 


= 5 . 5 OO = 2500 (flat plate) 


crit 

This critical Reynolds number must be compared with 


(13-4) 


U max r 


at the 


transition point for the pipe. Due to the parabolic lam i nar velocity 
distribution in the pipe = 2 u, and because r = d, then, for 

the pipe, i ^ HT r/v = ud/v. According to equation ( 13 . 1 ), = 2300 

for the pipe. Thus the comparable critical Reynolds numbers for pipe and 
plate show rather good agreement. _ 


CHAPTER XIY. OLDER THEORIES 

The first efforts toward theoretical calculation of the turbulent 
flows go back to Reynolds. One distinguishes in the theory of turbulence 
two main problems: 

1. The flow laws of the developed turbulent flow: 

The space and time velocity variations affect the time average of 
the velocity; they abt like an additional internal friction. The problem 
Is to calculate the local distribution of the time average of the velocity 
components, and thus to gain further information concerning, for instance, 
the friction drag. 
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2. Origin of turbulence: 

One Investigates under what conditions a small disturbance , 
superposed on a laminar .flow, increases with t ime . According to whether 
or not the disturbance increases with time, the l amin ar flow is called 
unstable or stable. The investigation in question is therefore a stability 
investigation, made to clarify theoretically the laminar— turbulent transi- 
tion. These investigations aim particularly at the theoretical calculation 
of the critical Reynolds number. They are, in general, mathematically 
rather complicated. 

The first problem, since it is the more important one for general 
flow problems will be our main concern. The second will be discussed 
briefly at the end of the lecture series. 

As to the first problem, that of calculation of the developed 
turbulent flow, one may remark quite generally that a comprehensive 
theoretical treatment, as exists for laminar flow, is not yet possible. 

The present theory of the developed turbulent flow must be denoted as 
semi— empirical. It obtains its foundations to a great extent from 
experiment, works largely with the laws of mechanical similarity, and 
always contains several or at least one empirical constant. Nevertheless 
the theory has contributed much toward correlating the voluminous experi- 
mental data and also has yielded more than one new concept. 

For the numerical treatment one divides the turbulent flow, unsteady 
in space and time, into mean values and fluctuation quantities. The mean 
value may a priori be formed with respect to either space or time. We 
prefer, however, the time average at a fixed point in space, and form 
such mean values of the velocity, pressure, shearing stress, etc. In 
forming the mean values one must not neglect to take them over a suffi- 
ciently long ti^pe interval T so that the mean value will be independent 
of T. Let the velocity vector with its three mutually perpendicular 
components be 


= iy + + kw (14.1) 

For a turbulent flow the velocity components are therefore functions of 
the three space coordinates and the time: 

u = u (x, y, z, t) 

v = v (x, y, z, t) (14.2) 

w = w (x, y, z, t) 

The time average for the component u, for instance, is formed as follows: 

1 / ,t ° 4T 

u(x, y, z) = - / udt (14.3) 

K 

■^Throughout the text, underscored letters are used in place of 
corresponding German script letters used in the original text. 
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If u, 7j w are independent of t Q and T, the motion is called steady 
on the average , or quasi steady. A steady turbulent flow, in the sense 
that the velocity at a point fixed in space is perfectly constant, does 
not exist. The velocity fluctuations are then defined by the equations 


u = u + u* 

V = V + V* i 

w = w + w’ 


(I*.*) 


and in the same way for the pressure: 


P = P + P’ 


(14-5) 


The time average of the 'fluctuation quantities equals zero, according to 
definition, as the following consideration will show immediately: 


Thus: 



u’dt 



n t o +T 


,, u 
u dt - - 


Jt 


dt = u — u = 0 


(14.6) 


u» 


V* 


= w 1 = 0 


(14.7) 


The Additional "Apparent" Turbulent Stresses 

i 

As a result of the velocity fluctuations additional stresses 
( = apparent friction) originate in the turbulent flow. This is readily 
illustrated for instance by the case of the simple shearing 'flow u =* u(y) 
(fig. 76). Here T = O’, however, a fluctuation velocity v* in the 
transverse direction is present. The latter causes a momentum transfer 
between the adjoining layers across the main flow. This momentum transfer 
acts like an additional shearing stress t. Whereas in laminar flow the 
friction is brought about by the molecular momentum exchange, the turbulent 
exchange of mome ntum is a macroscopic motion of, mostly, much stronger effect. 

The equations of motion of the turbulent flow, with this turbulent 
apparent friction taken into consideration, can be obtained from the 
Navier— Stokes differential equations by substituting equation (lk.4) into 
the latter and then forming the time averages in the Navier— Stokes 
differential equations. To that purpose the Navier— Stokes differential 
equations (3.I0) are written in the form: 
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8u C$(u 2 ) d(uv) d(uw) I 

[St + "^" + _ §r + — 




+ 






p 


f 8v 8(uv) 8 v 2 8(tw)1 

[8t 8x 8y 8z J 


8p 

8y 



+ 




p 


8w 8(wu) 
8t 8x 


d(wv) + dw 2 
8y 8z 


= - .§£ + n + &L + 

h * W 8/ 



> (1^.8) 


8u 

8x 


+ 




/ 


By introducing equations (14.4) and ( 14 . 5 ) and forming the time averages 
one first obtains from the continuity equation: 


and thus also: 





(14.9) 


8u* + 8v* + 8w* 
8x 8y 8z 


( 14 . 10 ) 


By introduction of equations (14.4) into the left side of equation (14.8) 
one obtains expressions as, for instance. 


u 


2 2 2 

= (u + u* ) = u + 2 u u* + u* etc. 


In the subsequent formation of the time average the squared terms in 
the barred quantities remain unchanged since they are already constant 
with respect to time. The mixed terms, as for instance u u*, . . . 
and also the termB that are linear in the fluctuation quantities are 
eliminated in forming the average because of equation (14.7). Howe ver, 

2 

the terms that are quadratic in the fluctuation quantities as u* , 
u*v f , . . . remain. Thus one obtains from the equation system 
(equation (l4.8)), after forming the time average the following system 
of equations: 
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w^ 

3z 







— 8w 
+ V 8y 






mAu 


l$u* 2 

s - 1 

+ Bu’v* 

, Su’w’ 

8x 


^ Sx 

'Sy 

8z 

_ + 

pAv 

-p \ 

8u' v* 

. + dv* 2 

! du’w* 

dy 

1 

5x 

Sy 

8z 

_ + 

pAv 

_r\ i 

Idu’v* 

, 8u'w* 

. + ^w* 2 

8z 

P 

8x 

Sy 

8z 


(14.11) 


The left side now formally agrees with the Ravier-Stokes differential 
equations for steady flow if one writes instead of u, v, w the time 
averages of these quantities. On the right side additional terms which 
arise from the fluctuations have heen added to the pressure and friction 
terms. 


Remembering that in deriving the Navier-Stokes differential equations 
one could write the resultant surface force per unit volume by means of 
the components of a stress tensor according to equation (3.7) in the form 


H = i 


d g x + ^ T xy + 5t xz 
^ 8x §y 8z j 



(14.12) 


one recognizes by comparison with (equation (l4.11)) that one may introduce 
for the quantities added by the fluctuation motion a symmetrical stress 
tensor in the following manner: 


ff x 

.2 

= — p u' 

T 

xy 

= — pu * V 1 

T 

xz 

= — pu’w’ 


; 


T 

xy 

= — pu*v* 

°7 

II 

1 

TO 

tl 

rol 

T 

yz 

= — pv*w* 

. 

1 


(14.13) 

T 

xz 

= — pu'w* 

T 

yz 

= — pv’w’ 

ff z 

= — pw* 2 





-/ 


One has therefore, for the mean values of the quasi s beady flow, the 
following equations of motion: 
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p 




p 





dp 

— + pAu + 

OX 


dcx„ dr 


X 

dx 


xy 


3t 


dy 


xz 

dz 


s- St 

_ k + ^7 + -S 

dy 


dcr dT 


+ — - + 

dx dy dz 


yz 


> (14.14) 



v§£ + 
dy dz / 


. dP 
• dz 


+ 


pAw + 


dT 


xz 


dx 


dT 


_yz 


dy 



y 


The continuity equation (equation ( 14 ■ 9 ) ) also enters. The boundary 
conditions are the same as for laminar flow; adhering of the fluid to the 
wall, that is, on the solid walls all velocity components eq ual zero. 
According to equation (14.14) the mean values of the turbulent flow obey 
the same equations of motion as the velocity components of a 1 nmlrmr flow, 
with the friction forces, however, increased by the apparent stresses of 
the turbulent fluctuation motion. But since the fluctuation velocities 
} • • « and particularly their space distribution are unknown, 
equations (14.14) and (14.13) are, at first, rather useless for the 
calculation of a turbulent flow. 


Only when one will have succeeded in expressing the fluctuation 

p — - 

quantities u’ , u’v’, ... in a suitable manner by the time averages 

u, v, . . ., will it be possible^to use equations (14.14) to calculate, 
in particular, the mean values u, "v, "w. 

A first expression of this kind which brought however little success 
was originated by Boussinesq (reference 64). He introduced, aside from 
the ordinary viscosity coefficient, a new viscosity coefficient of the 
apparent turbulent friction. In analogy to the stress tensor for 
laminar flow which Is, according to equation (3.13): 


T xy T xz 


T xy °y T yz 


T xz T yz °z 


p 0 0 
0 p 0 
0 0 p 


+ H 


du du du 


du dv dw 

dx dy dz 


dx dx dx 

dv dv dv 

-1- 1 1 

du dv dw 

dx dy dz 

-h |J. 

dy dy dy 

dw dw dw 


du dv dw 

dx dy dz 


dz dz dz 


S 55 (14.15) 


Boussinesq puts for the apparent turbulent friction: 


■ 2p 6 S V ” p 6 


( du , dv\ 
\dy d z y 


(14.16) 


Then there corresponds to the laminar viscosity coefficient p the mixing 
factor p g; 
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p ~ p e or V~e 


The kinematic viscosity of the turbulent flow (apparent friction) is 
usually very much larger than that for the ordinary laminar friction. 
(Hundred— or thounsandfold or more). In general, one may therefore 
altogether neglect the ordinary viscosity terms pAu, ... in 
equation (lk.l4). Only at the solid walls where due to the no— slip 
condition - — 


u=v=w=0 as well as u* = v* = w’ = 0 


the apparent turbulent friction disappears, does the laminar friction 
again become dominant. Thus there exists in every turbulent friction 
layer in the Immediate neighborhood of the wall a very narrow zone where 
the flow is laminar. The thickness of this laminar sublayer is only a 
small fraction of the turbulent boundary layer thickness. 

One can easily understand from the example of the simple s hear flow 
according to figure 7 6 that in a turbulent flow the mean value u’v* is 
different from zero. For this case, a correlation exists between the 
fluctuation velocities u* and v* in the following manner: The 

particles with negative v* . have "mostly** a_positive u', since they 
come from a region of larger mean velocity u.* The parts with positive 
v', on the other hand, have "mostly" a negative u’, because they come 
from a region with smaller "u, and retain in the transverse motion 
approximately the x— momentum of the layer from which t hey come. Thus, 
'"mostly" u’v* < 0 and, therefore, the time average u* v* < 0. Therefore, 
the shearing stress Is this flow Is: 


t = — p u’v' > 0 
xy p 


In measuring turbulent flows one usually measures only the mean 
values u, v, . . • . since only they are of practical interest. However, 
in order to obtain deeper insight into the mechanism of the turbulent 
flow, the fluctuation quantities have recently been measured and also 
their mean squares and products: 

Vu « 2 \Zv ’ 2 \/u*v« 


According to measurements by Reichardt (reference 6 5) in a rectangular 
tunnel (width 1 millimeter, height 2k centimeters) the maximum value 


* "Mo ably" is to indicate that particles with different signs,, 
thouph not excluded, are in the minority. 
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of V u * 

equals O.Op u WflT . Both maxima lie in the neighborhood of the wall. 

One may say, therefore, that in this case the turbulence is strongest 
near the wall. 

In a flow that is homogeneous (wind tunnel), turbulent fluctuations 
are also always present to a varying degree. They determine the so-called 
degree of turbulence of a wind tunnel. Since the measurement of the 
fluctuation quantities is rather difficult (hot-wire method), a more 
convenient measuring method has been chosen, for the present, for determin- 
ation of the degree of turbulence of a wind tunnel: namely, the determin- 

ation of the critical Reynolds number for the sphere from force measure- 
ments or pressure distribution measurements. One defines aB critical 
Reynolds number the one where the drag coefficient c^. = 0.3. It 

becomes clear that a unique connection exists between the critical Reynolds 
number and the turbulent fluctuation velocity in the sense that the 
critical Reynolds number of the sphere is the lower, the higher the 
turbulent fluctuation velocity. According to American measurements 
(reference 97) the connection between the longitudinal fluctuation and 
the measured critical Reynolds number for the sphere is as shown in the 
following table: 


, for instance, equals 0.13 u ma ^, the maximum value of V v 


\t- 


& 

0.004 

0.0075 

0.012 

0.017 

0.026 

Re ,, 10-5 
crit 

2.8 

2.4 

2.0 

1.6 

1.2 


In addition to the apparent increase of viscosity, the turbulent 
fluctuation motion has other effects: It tends to even out any tempera- 

ture differences or variations in concentration existing in a flow. 

The diffusion of heat, for instance, is much larger than for laminar flow, 
because of the exchange motions which are much stronger in turbulent 
flow. A close connnection therefore exists, for instance, between the 
laws of flow and of heat transfer from a heated body to the fluid flowing 

fcy* 


Ninth Lecture (February 2, 19*1-2) 
CHAPTER XV. MORE RECENT THEORIES; MIXING LENGTH 


In order to make possible a quantitative calculation of turbulent 
flows, it is necessary to transform the expressions for the apparent 
turbulent stresses (equation (14.13)) in such a manner that they no longer 
contain the unkn own fluctuation velocities but contain the components of 
the mean velocities. Consider, for that purpose, a particularly simple 
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flow, namely a plane flow which hae the same direction everywhere and a 
velocity varying only on the different stream lines. The main— flow 
direction coincides with the x— direction; then: 


u = u(y) v = 0 w = 0 


(15.1) 


Of the shearing stresses., only the component t = t is present, for 

which from equation (14.13) as well as fromBoussinesq’s equation, 
equation (14.15), there results: 


T 


p U*V* = P € — 

d7 


(15.2) 


This formula shows that | t |/ p equals the square of a. velocity. 
putSj therefore, for use in later calculations. 



One 


(15.3) 


and denotes v* as shearing stress velocity. Thus this shearing stress 
velocity is a measure of the momentum transfer "by the turbulent 
fluctuation motion. 

According to Prandtl (reference 66), one may picture the turbulent 
flow mechanism, particularly turbulent mixing, in the following simplified 
manner: Fluid particles, each possessing a particular motion, originate 

in the turbulent flow; they move for a certain distance as coherent masses 
maintaining their velocity (momentum) . One now assumes that such a fluid 
particle which originates in the layer (y^ — l) and has the velocity 
ufy^ — l) moves a distance l = mixing length normal to the flow (fig. 77). 
If this fluid particle maintains its original velocity in the x-direction 
it will have, in its new location yj_, a smaller velocity than its new 
surroundings, the velocity difference being 

-u| = u(y 1 ) - u(y 1 - Z) with v* > 0 


Likewise a fluid particle coming from the layer (y^ + Z) to y^_ has 
at the new location a greater velocity than the surroundings there; the 
difference is 


n 
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u 2 = U-(72_ + 1) - with v* < 0 


u| and give the turbulent velocity fluctuation in the layer y . 

One obtains for the mean value of this velocity fluctuation 


|U'| = | (|u* I + fujl) 



( 15 . 4 ) 


From this equation one obtains the following physical interpretation for 
the mixing length 2 : 

The mixing length signifies the distance in the transverse direction 
which a fluid particle must travel at the mean velocity of its original 
layer so that the difference between its velocity and the velocity of the 
new location equals the mean velocity fluctuation of the turbulent flow. 

It is left open whether the fluid particles in their transverse motion 
fully mainfain the velocity of their original layer, or whether they have 
partly assumed the velocity of the traverse layer and then travelled 
larger distances in the transverse direction. The Prandtl ml Ting length 
which is thereby introduced has a certain analogy to the mean free path 
of the kinetic theory of gases, with, however, the difference that there 
one deals with microscopic motions of the molecules, here with macroscopic 
motions of larger fluid particles. 

One may picture the origin of the transverse fluctuation velocity v’ 
In the following way: 

Two fluid particles flowing from the layers (y^^ + l) and (y^ — 2) 

meet in the layer jj_ in such a manner that one lies behind the other: 
the faster (y-j_ + Z) behind the slower (y^ — 2). They then collide 

with the velocity 2u* and giveway laterally. Thereby originates the 
transverse velocity v' 1 , directed away from the layer y^ to both sides. 

If, conversely, the slower of the two particles is behind the faster, 
they withdraw from each other with the velocity 2u* . In this case the 
space formed between them is filled up out of the surroundings. Thus 
originates a transverse velocity v* directed toward the layer y^_. 

One concludes from this consideration that v* and u’ are -of the same 
order of magnitude and puts 


j v 1 1 = number 


|u’| = number Z — 

<3y 


( 15 - 5 ) 


Tn order to express the shearing stress according to equation ( 15 . 2 ) one 
has to consider the mean value u*v* more closely. The following 
conclusions can be drawn from the. previous considerations. 
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The particles arriving In the layer y-^ with positive v* (from 

"below, fig. 77) have "mostly" a negative u* so that u’v’ is negative. 
For the particles arriving with negative v’, u* is "mostly" positive, 
so that u’v’ is again negative. "Mostly" signifies that particles with 
different sign a re n ot wholly excluded, "but are strongly outnumbered. 

The mean value u’v* is therefore different from zero and negative. Thus 
one puts 


u’v’ = — k |u’| I v ’ | 

with k / Oj 0 <k < 1. The numerical factor k, 
correlation coefficient is not known more closely, 
equation (15,5) and (15.^) one now obtains 

u’v’ = — number ] 

w ; 

the "number" in this equation being different from the one in 
equation (15.5). If one includes the "number" in the unknown mixing 
length, one can also write 


( 15 . 6 ) 


also called the 
According to 


(15-7) 


u’v’ = - 



(15.8) 


and thus finally obtains for the turbulent shearing stress according to 
equation (15.2) 


T = 



Considering that the sign of t also must change with the sign of ^ 
It is more correct to write 


cu 

r-5 

a 

It 

t- 

du 

du 

Prandtl ’ s 



dy 

formula 


(15.9) 


This is the famous Prandtl mixing lfength formula which has been very 
successful for the calculation of turbulent flowB. 


If one compares this formula (equation (15.9)) with the equations 
of Boussinesq where one had put t = e (e = mixing factor = turbulent 
analogue of the laminar viscosity a), erne has for the mixing factor 
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e 



du. 

ay 


(15.10) 


The turbulent mixing factor e is in most cases larger than the laminar 
viscosity 4 by several powers of ten. Moreover the mixing factor e 
is dependent on the velocity and on the location and tends toward zero 
near a wall, because there the mixing length goes toward zero. 

If one compares Prandtl’s formula equation (I5.9) with Boussinesq * b 
equation ( 15 . 2 ) one could perhaps think at first that not much has been 
gained, since the unknown quantity e ( = apparent viscosity) has been 
replaced by the new unkn own 2 = mixing length. Nevertheless Prandtl’s 
formula is considerably better than the 'old formula for the following 
reason: It is known from tests that the drag for turbulent flow is 

proportional to the square of the velocity. According to equation (15.9) 
one obtains this square law for drag by assuming the mixing length to be 
independent of the velocity, that is, by assuming the mixing length to 
be purely a function of position. It is considerably easier to make a 
plausible assumption for the length 2 = mixing length than for the 
apparent turbulent viscosity e , and therein lies the considerable 
superiority of Prandtl’s formula equation ( 15 . 9 ) over Boussinesq’s 
equation (15.2). 

In many cases the length 2 can be brought Into a simple relation 
to the characteristic lengths of the respective flows. For the flow along 
a smooth wall 2 must, at the wall itBelf, equal zero, since all trans- 
verse motions are prevented at the wall. For the flow along a rough wall, 
however, the limiting value of 2 at the wall equals a length of the 
order of magnitude of the height of the roughness. 

It would be very useful to have a formula permitting the determina- 
tion of the dependence of the mixing length on the position for any 
arbitrary flow. Such an attempt has been made by v. Karman (reference 68). 
v. Karmen makes the assumption that the inner mechanism of the turbulent 
flow is such that the motion at various points differs only with respect 
to t ime — and length-scale, but is otherwise similar (similarity hypothesis). 
Instead of the units of time and length one may select those of velocity 
and length. The velocity unit that 
is the shearing stress velocity v* 
corresponding unit of length Is the 

Tn order to find the quantity 
u(y), Karman applies the Taylor 
borhood of the point y 1 - 

U(y) = u( 7l ) + (y - y x ) (^) x 

* In the following, the bar over the mean velocity will be omitted, 
for simplification. 


is important for the turbulent motion 
according to equation ( 15 » 3 )» The 
m ix ing length 2 . 

2 from the data of the basic flow 
development* for u(y) In the neigh- 


+ Fr( y - sr i )2 ( 4 ¥) + (15,l:L) 
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The length 2 cannot depend on the velocity u(y^), since according 

to Newton’s principle of relativity the addition of a constant velocity 
has no influence on the course of motion. Thus 



*s 


\ 


and the higher derivatives remain as characteristic data of the hasic 
flow. The simplest length to he formed from it is 



v. Karman puts therefore 


du 

dy 



( 15 . 12 ) 


According to this formula 2 is not dependent on the amount of velocity 
hut only on the velocity distribution. Thus 2 is a pure position 
function as required above. In equation (15.12) is an empirical 

constant which must he determined from the experiment. To arrive further 
at the turbulent shearing stress, v. Karman also maintains Prandtl’s 
equation (I 5 . 9 ). 


In generalizing equation (15.9) one obtains, according to Prandtl, 
the complete expression for the turbulent stress tensor of a plane flow 
in the form 



(15.13) 




The common factor on the right side signifies the turbulent mixing 
factor according to equation (I 5 .IO). 
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Tenth. Ledture (February 9 , 19^2) 


Flow Along a Smooth Wall 


We will immediately make a first application of Prandtl's formula 
(equation (I5.9)) for the flow along a smooth wall. The normal distance 
from the wall is denoted as y. Let the wall coincide with the x— axis. 
For the Telocity distribution then, u = u(y). For this case one sets 
the mixing length in the neighborhood of the wall proportional to the 
distance from the wall 


2 = kj (15. 1*0 

the constant k must be determined from the experiment. Moreover one 
makes the assumption that the shearing stress t is constant in the 
entire flow region; then the shearing stress velocity v* according to 
equation (15*3) also is constant. If one further neglects the laminar 
friction, one obtains from equations (15* 2), > (15*9 ) » (I5*b^) 


v* 


2 = 



or 


du _ v* 

dy rcy 


and by integration 


v* 

u = — Zn y + constant 


(15.15) 


In determining the constant of integration one must pay attention to the 
fact that the turbulent law equation (15-9) does not apply right up to 
the wall but that very near to the wall an extremely thin laminar layer 
is present. From the laminar viscosity p and the turbulent shearing 
stress velocity v* one can form the length v/v*. The constant of 
integration in equation (15.15) is determined from the condition that 
u = 0 for y = y Q . Thus there results, according to equation (15.15) 

u j (In y - In y°) (15-16) 


The as yet unknown distance from the wall y D is set proportional to 
the length v/v*, thus 
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7o = P ^ (15.17) 

where j3 signifies a dimensionless constant. Thus one finally obtains 
for the velocity distribution at the smooth wall 

u = TT T* ~ 2n p) (15.18) 

that is, a logarithmic velocity distribution law. It contains two 
empirical constants n and j3. According to measurements k = 0.4. 

From equation (15.18) one can see that the dimensionless velocity u/v* = cp 
can be represented as a function of the dimensionless distance from the 
wall T) = v*y/v. The latter is a sort of Reynolds number, formed with the 
distance from the wall y and the shearing stress velocity v*. Thus one 
obtains for larger Reynolds numbers from equation (15.18) the following 
universal velocity distribution law 


cp(tj) = A 2n tj + B 


(15.19) 


with A = l/n - 2.5. For smaller Reynolds numbers, where the laminar 
friction also has a certain influence, tests gave the velocity distribution 
law 


or 


9(tj) = C i\ 


n 


(15.20) 



(15.20a) 


with the exponent n equalling about l/7« These universal velocity dis- 
tribution laws according to measurements for pipe flow are given in 
figure 78. They will be discussed in more detail in the following chapter. 


CBAPTSR XVI. PIPE FLOW 
a. The Smooth Pipe 


Among the various turbulent flows of practical importance, pipe 
flow was investigated with particular thoroughness because of its great 
practical importance. We shall therefore consider the pipe flow first. 

It will be noted at this point that the flow laws of the pipe flow may 
be applied to other cases, as for instance the plane plate in longitudinal 
flow. Consider a straight pipe of circular cross section and with a smooth 
wall. Let y be the radial coordinate measured from the pipe axis. The 
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"balance of forces "between the shearing stress t and. the press'ore 
drop p^ — pg on a piece of pipe of length L yields as "before for the 
laminar flow according to equation (2 ►la), the relation: 


T _ p i ~ p 2 z 

L 2 


(16.1) 


This for mul a applies equally to laminar and turbulent flow. In it t 
now signifies the sum of the laminar shearing stresses and of the apparent 
turbulent shearing stresses. Over a cross section, t is proportional 
to y. The shearing stress at the wall t 0 may be determined 
experimentally by measurement of the pressure drop: 


T 


O 


Pi ~ *2 r 
L 2 


(16.2) 


For the turbulent flow the connection between pressure drop and flow 
■volume Q = jrr 2 u must be obtained from tests.* In the literature there 
exists a very great number of pipe resistance formulas. Only those serve our 
purpose which satisfy Eeynolds* law of similarity. One of them is the 
formula of H. Blasius (reference 69), set up particularly carefully, 
which is valid for a smooth wall and for Eeynolds numbers 
Ee = ud/v ^ 100 000. 

If one introduces, as before in equation (2.6), the dimensionless 
pipe resistance coefficient X by the equation 



(16.3) 


X 


is, according to Blasius: 

\ = 


0.3l6h 



(1 6.k) 


Comparing equations (16.2) and (16.3) one finds: 


X 

S 


T . = T p u 


and therefore according to equation (1 6.4): 


(16.5) 


T Q = 0.03955 p jA V lA d (16.6) 


*In the following, u is, for the pipe flow, the mean flow velocity 
at the cross section, as distinguished from the time average in the previous 
sections. 
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If one introduces, moreover, instead of the diameter d the radius r, 
the numerical factor in this linear equation must he divided by 

2^A = 1.19* Thus r 0 becomes: 


T o = O.O3325 p u T A v 1 / 1 *- r^A = p v* 2 (16.7) 


where the shearing stress velocity is defined by the wall shearing stress: 


v* = 



To 


(16.8) 


p 7 /lj. 

If one finally factors the quantity v* in equation (16.7) into v* ' x 


1/4 


one obtains: 




Y/7 


(i6r 9 ) 


This equation is very similar to equation ( 15.20a); however, the mean 
velocity now takes the place of the local velocity and the pipe radius 
takes the place of the distance from the wall. One passes first from 
the mean velocity to the maximum velocity u^ ; based on measurements of 

Nikuradse (reference 70) u = 0.8 u^, and therewith follows from 

equati on ( 16 . 9 ) : 


If this formula is assumed to be valid for any distance from the wall, 
one obtains: 


* • »■» (¥f 


or 


• V / 
<p = 8.74 


(16.10) 


(16.11) 


This is the so-called 1/7— power law for the velocity distribution; its 
form was already given in equation (15.20a). The coefficients n and 
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still unknown there, have now "been determined on the "basis of the 
resistance law of the pipe flow. Figure 78 shows, according to measure- 
ments of Nilruradse (reference 70} that this law is well satisfied in the 
range of Reynolds numbers up to 100,000. Naturally this law of velocity- 
distribution can apply only to the region of Reynolds numbers where the 
pipe resistance law given by equation (l6.4) is valid, since it was 
derived from this law. 

For purposes of later calculations we shall derive from 
equation (16.10) the shearing stress velocity v*. One obtains: 


7/8/ v x 1/8 

V* = 0.150 u (Jrj 


(16.12) 


with 8.7h 7//8 = 6.65 and 


6.65 


= 0.150. From equation (16.12) follows: 


T o = P T 


, 2 = 0.0225Pu 7//if ' 


(16.13) 


This relation will be needed later. 

Comparing measured velocity distributions with equation (16.10) 
one can state that outside of the range of validity of equation (16.10), 
namely for Re >100,000, a better approximation 1 b obtained by the 
power l/8, 1/9, or l/lO instead of l/7» The measurements concerning 
the pipe resistance (fig. 8l) show an upward deviation from the formula 
of Blasius for larger Reynolds numbers. 

The logarithmic velocity distribution law, equation (15.19), derived 
in the previous chapter has been verified by Nilruradse (reference JO) on 
the basis of his measurements for the smooth pipe. For this purpose from 
the measured pressure drop for each velocity profile one first deter min es 
the wall shearing stress according to equation (16.2) and from tha t, 
according to equation (16.8), the shearing stress velocity v* = \J r Q /p- 

Then the dimensionless velocity cp = u/v # can be plotted against the 
dimensionless distance from the wall q = yv*/v. The measurements of 
Nilruradse in a very large range of Reynolds numbers. Re = 4 X 10^ up to 
3240 X 1C>3, lie very accurately on a straight line if one plots cp against 
log Tj (fig. 78). The straight line has the equation; 


cp = 2.5 In t) + 5.5 


(16.14) 


This gives, by comparison with equation (15.18), the following numerical 
values for the coefficients K and p 
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k = o.koo 3 = o.m 


( 16 . 15 ) 


Mixing Length 


From the measured Telocity distribution and the measured pressure 
drop the distribution of the mixing length over the pipe cross section 
can be determined according to equations ( 16 . 2 ), ( 16 . 1 ), and ( 15 . 9 ) • 

T = t q ~ (y = distance from the pipe axis). This determination of the 

mixing length from the measurements in the pipe was made by Nikuradse 
(reference 70). For large Reynolds number s, where the influence of 
viscosity is negligible, one obtains a distribution of the mixing length 
if r over y/r which is independent of the Re-number (fig. 79)* The 
following interpolation formula can be given for this distribution: 

1 = o.lk - 0.08 f l - £) - o.o 6 [i - 2y (16.16) 

In this equation y signifies the distance from the wall. The develop- 
ment of equation (l6.l6) for small y/r (neighborhood of the wall) gives 

y 2 

l m O.ky - O.kk ~ + . > . (l6.l6a) 


In the neighborhood of the wall the mixing length is, therefore, propor- 
tional to the distance from the wall. Equation (l6.l6) for the distri- 
bution of the mixing length applies not only to the smooth pipe, but, 
according to the measurements of Nikuradse (reference 71) also to the 
rough pipe, as can be seen from figure 79* From this fact one can derive 
in a very simple manner a universal form for the law of velocity distri- 
bution, valid for the smooth as well as for the rough pipe. One puts 

for the mixing length distribution: l = kj f^with f ^)— > 1 for 2. — > 0. 

Furthermore follows from the linear distribution of the shearing stress 
over the cross section: 


T 



(y = distance from the wall) 


together with equation ( 15 . 9 ) 


du _ 1 .ff~ 
dy 2 Y p 



K yf(y/r) 


( 16 . 17 ) 


and hence by integration: 
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(16.18) 


the lower l imi t of intergration y Q where the velocity equals zero 1 b, 
according to the considerations of the previous section, proportional 


to 


v/v*; thus: 



From equation (l6.l8) follows: 



V 1 - y/r d. y/r 



(16.18a) 


and therefore, from equations (l6.l8) and (l6.l8a): 


u^ “ u = T * F 


(16.19) 


This law, with the same function F(y/r), applies equally to smooth and 
rough pipes. Xt states that the curves of the velocity distribution over 
the pipe cross section for all Reynolds numbers and all roughnesses can 
be made congruent by shifting along the velocity axis, if one plots 
u mpT — u/v» against y/r (fig. 80). The explicit expression for the 

function F(y/r) is obtained immediately from equation (l6.1h), according 
to which 


^ - u =t 2.5v* In 21 = 5.75 -7* log (16.20) 

y y 


Universal Resistance Law 


.A 




According to their derivation the velocity-distribution— law 
(equations (16.19) and (16.20)) are to be regarded as valid for arbitrary 
Reynolds number since the laminar viscosity was neglected as compared 
with the turbulent viscosity. We shall now derive from the velocity- 
distribution— law equation (16.20) a resistance law which, in contrast to 
Blasius*, applies up to Reynolds numbers of arbitrary magnitude. 
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From equation (16.20) one may determine by integration over the 
cross section the mean flow velocity u. One finds: 


u ~ u max 3 * 75 v* 


(16.21) 


The test results of Nikuradse -(reference 70) gave a number slightly 
different from 3.75, namely: 


* = Vi “ k -° 7 T * 


(16.22) 


According to equation (16.5): 


X 



(16.23) 


From the universal velocity distribution law of the smooth pipe 
equation (l6.1k) follows: 


^ rnn.Y 


|2.5 2 n 


rv* 

V 


+ 



and hence the connection with equation (16. 2l): 


u 


{■ 


2.5 2n 


rv* 

V 


+ 



(16.2k) 


The Reynolds number enters into the calculation by means of the identity: 


1 ud v* _ ud Vx 
V " 2 V T " V }±\f~2 


Thus results from equations (16.23) and (16.2k): 


X = 


8 


{"2.5 in (¥ vj) - 2.5 In k \[2 + 1.75] 
' * 


\ 


2.035 log(^ Vx)-0.9lj 
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<^*1 


-r 


or: 


i = 2.035 log 

ft 



(16.25) 


Accordingly a straight line must result for the resistance law of the smooth 

pipe, if one plots 1 / against log^pVT^. This is very well confirmed 

hy Nikuradse’s measurement (fig. 8l) . • The numerical values according to the 
measurements differ only slightly from those of this theoretical derivation. 
From. Nikuradse’s measurements was found: 



(16.26) 


Universal Resi stance Law For Smooth Pipes 

This is the final resistance law for smooth pipes. On the "basis of 
its derivation it may he extrapolated up to Reynolds numbers of arbitrary 
magnitude. Thus measurements for larger Reynolds numbers than those of 
Nikuradse's tests are not required. Up to Re = 100,000 this universal 
resistance law is in good agreement with the Blasius law according to 
equation (l6.li-). For higher Reynolds numbers the Blasius law deviates 
considerably from the measurements (fig. 8l). 

Concerning the determination of X, from equation (16.26), where it 
appears on both sides, it should be added that it can be easily obtained 
by successive approximation. 


A 




Eleventh Lecture (February 16, 19^2) 


b. The Rough Pipe 


The characteristic parameter for the flow along a rough wall is the 
ratio of grain size k of the roughness to the boundary layer thickness, 
particularly to the thickness of the laminar sublayer 5^ which is 


always present within the turbulent friction layer in the immediate 
neighborhood of the wall. The thickness of the laminar sublayer is 


6^ = number — . The .effectiveness of roughness of a certain grain size 
depends, therefore, on the dimensionless roughness coefficient k/8 ? ~ 


kv* 


In the experimental investigations of the resistance of turbulent flows 
over rough walls, the rough pipe has been studied very thoroughly since it 
is of great practical importance. Besides depending on the Reynolds number, 
the resistance of a rough pipe is a function of the relative roughness r/k. 
One distinguishes for the resistance law of a rough pipe three regions. 
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Tlie subsequently given "boundaries of these regions are valid for sand 
roughness k 0 like those investigated by Nikuradse (reference 7l). 


1. Hydraulically smooth: Th9 grain size of the roughness is so 
nrnn.1 1 that all roughnesses lie within the laminar sublayer. In this 
case the roughness has no drafe increasing effect. This case exists for 
small Reynolds numbers and for values of the characteristic roughness 

< T ** a < 


number: 


0 ^ _U3 5. 

v 


2. Fully developed roughness flow; The grain size of the roughness 
is so large that all roughnesses project from the laminer sub-layer. The 
friction drag then consists predominantly of the form drag of the single 
roughness elements. A purely square drag law applies. _ For the pipe the 
resistance coefficient X is then independent of Re and only dependent 
on the relative roughness k/r. This law exists for^vgry large Reynolds 

numbers. For sand roughness tills law applies for: — ' S ^ TO. 


3. Intermediate region: Only a fraction of the roughness ‘elements 
project from the laminar sublayer. The drag coefficient depends on r/k 
as well as on Re. This law exists for medium Reynoldsjaumbers and, for 


the sand roughness, for: 


5 — 


I^< 70. 


The dependence of the pipe resistance coefficient on the Reynolds 
number and on the relative roughness according to the measurements of 
Nikuradse (reference 71) can be seen from figure 82 as well as, in 
particular, the three laws just given. 


The velocity distribution on a rough wall is given, basically, in 
the same way as for the smooth wall by equation (15.16). One has only 
to substitute for the constant of integration y Q another value: y D 

proportional to the roughness grain size. One puts for sand roughness 
y Q = 7 k 0 and hence obtains from equation (15.16) 



r 





2n 7 ► 


(16.27) 


A 




The constant 7 is, moreover, a function of the roughness form and the 
roughness distribution. Comparison with experiments of Nikuradse 
(reference 71) on pipss roughened artificially by sdnd yields for the 
velocity distribution the general formula: 

u = v # ^2.5 2n ~ + B 




(16.28) 



NACA TM No. 1218 


27 


7 


the constant B being different in each domain described above; it 
depends on v # k s /v . 

For the fully developed roughness flow B = 8.5, thus: 

u “ In p- + 8.5^ (fully rough) (16.29) 


whence follows: 


and: 


W = T *( 2 '5 * n ^ + 8*5) 

u max “ u = T * 2 * 5 - ^ | 


(16.30) 


(16.30a) 


in agreement with equation (16.20). Thus there applies also to rough 
pipes, as equation (16.21) did before to the smooth pipe, the relation: 


u = % ax ~ 3*75 v* 


(16.31) 


From here one can, by a calculation which is perfectly analogous to the 
previous one for the smooth pipe, easily arrive at the resistance law of 
the rough pipe for fully developed roughness flow. By insertion of 

according to equation (16.30) into equation (16.31) one obtains: 


or; 


or: 


u = v* 




(16.32) 




8 


(2.5 to X. + 4.75 J 


( 16 . 33 ) 


X = 


^2.0 log X + 1.68 J 


( 16 . 34 ) 
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This is the square resistance law of the fully developed "roughness flow. 
Comparison with the test results of Nikuradse (fig. 83) shows that one 
obtains better agreement if one changes the number 1.68 to 1.74. Thus 
the resistance law of the pipe flow for fully developed roughness is: 



log + 1.74 
*s 


(16.35) 


In the plots of log r/k S( (fig* 83) the test results fall 

very accurately on a straight line. 

For flow along a rough wall in the intermediate region the constant 
B in equation (16.28) is, moreover, a function of the roughness coeffi- 
cient v*k 0/ ^. For this case also one can derive the resistance law 

immediately from the velocity distribution. According to equation (16.28): 

B = — - 2.5 in -£• - - 2.5 In -£ (16.36) 

v^ k v„ k 


On the other hand, according to equations (16.31) and (16.23): 


iS.it 3-75 = 3.75 

T * T * v* 


(16.37) 


so that one obtains from equation (16.36): 

b(— ^r) = — - 2.5 to -2- - - 2.5 to p- + 3.75 (16.38) 

\ v ) v* k 0 “a 

One can, therefore, determine the constant B as a function of v*k^/v 

either from the velocity distribution or from the resistance law. The 
plot in figure 84 shows good agreement between the values determined by 
these two methods. At the same time the determination of the resistance 
law from the velocity distribution is confirmed. 

The formula for B includes the case of the smooth pipe. B is, 
according to equation (l6.l4). 
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B 




v*k 

2.5 2n + 5.5 

V 


(16.39) 


Thus a straight line results for B in the plot against log v*k 0 /v. 


Other Roughnesses 

Because of the great practical importance of the roughness— pr oh lem 
a few data concerning roughnesses other than the special sand roughness 
will he given. Wikuradse’s sand roughness may also he characterized hy 
the fact that the roughness density was at its maximum value, because 
the wall was covered with sand as densely as possible. For many practical 
roughnesses the roughness density is considerably smaller. In such cases 
the drag then depends, for one thing, on form and height of the roughness, 
and, moreover, on the roughness density. It is useful to classify any 
arbitrary roughness in the scale of a standard roughness. Nikuradse’s 
sand roughness suggests itself as roughness reference (roughness scale) 
because it was investigated for a very large range of Reynolds numbers 
and relative roughnesses. Classification with respect to the roughness 
scale is simplest for the region of fully developed roughness. According 
to what was said previously, for this region the velocity distribution 
is given by: 


£ = 5.75 l°& ^ + b b b b = 8.5 (16.40) 


and the resistance coefficient by: 


X - j ~ (16.41) 

f 2.0 log + 1.74 ) 


One now relates fco an arbitrary roughness k an equivalent sand roughness 
k B by the ratio 

k g = a k (l6.42) 

Where by equivalent sand roughness k B Is meant that grain size of sand 
roughness which has, according to equation (l6.4l) the same resistance 
as the given roughness k. 

Basic ally , of course, the equivalent sand rou ghn ess k B can be 
determined by a resistance measurement on the pipe. However, such measurements 
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for arbitrary roughnesses are difficult to perform. Measurements on *v 

arbitrary roughnesses in a tunnel with plane walls are more convenient. 

To this purpose an exchangeable wall of a tunnel with rectangular cross 
section is provided with the roughness to be investigated (fig. 85). 

From the measurement of the velocity distribution in such a tunnel with I"' 

a rough and a smooth longitudinal wall one obtains, for the logarithmic 
plot against the distance from the wall, a triangular velocity distribu- 
tion (compare fig. 85). From the logarithmic plot of the velocity distri- 
bution over the rough wall - - 


u = iij, log y + m^ 


(16.43) 


one obtains by comparison with the universal law according to 
equation (16.28) for the shearing stress velocity at the wall: 


T *r 


^r 

5-75 


(16.44) 


Further, one determines for the roughness to be investigated the constant 
B of the velocity distribution law, namely: 

B = — - 5.75 log | (16.45) 

T *r * 


By comparison of equation (16.45) with (16.40) one obtains for the 
equivalent sand roughness-; 

k_ 

5.75 log — = 8.5 -B (16.46) 

k 


In this way one may determine the drag for arbitrary roughnesses from a 
simple measurement in the roughness tunnel. This method may be also 
carried over to the case of the intermediate region. 


CHAPTER XVII. THE FRICTION DRAG OF THE FLAT PLATE 
IN LONGITUDINAL FLOW 


The turbulent friction drag of the plate in longitudinal flow is 
of very great practical importance, for instance as friction drag of wings, 
airplane fuselages, or ships. The exact measurement of the friction drag 
for the large Reynolds numbers of practice is extremely difficult. Thus 
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it is particularly important that one can, according to Prandtl 
(references 73 and 7^)s calculate the friction drag of surfaces from the 
results of pipe flow studies. This conversion from the pipe to the plate 
can he made for the smooth as well as for the rough plate. 


a. The Smooth Plate 


One assumes, for simplification, that the boundary layer on the 
plate is turbulent from the leading edge. Let the coordinate system be 
selected according to figure 86. The boundary layer thickness 5 (x) 
increases with the length of run x. Let b be the width of the plate. 
For the transition from pipe to plate the free stream velocity TT 0 of 

the plate corresponds to the maxi mum velocity in the pipe, and 

the boundary layer thickness 5 to the pipe radius r. 


One now makes the f undame ntal assumption that the same velocity 
distribution exists in the boundary layer on the plate as in the pipe. 
This is certainly not exactly correct since the velocity distribution in 
the pipe is influenced by a pressure drop, whereas on the plate the 
pressure gradient equals zero. However, slight differences in velocity 
distribution are insignificant since it is the momentum integral which 
is of fundamental importance for the drag. For the drag W(x) of one 
side of the plate of length x, according to equations (10. 1) and (10.2) 




V(x) = b 


a5(x) 


T n (x) dx = bp 


u (°o - u ) ay 


do 


(17.1) 


whence 


1 dW 
b dx 



(17.1a) 






The equation (17.1) can also be written in the form 

W(x) = T>pTT 0 2 8(i/ i (l - d | (17.2) 

For the velocity distribution in the boundary layer one now assumes 
the l/7— power law found for the pipe. Replacing u T11flT by U 0 and r 

by 5q one may write this law, according to equation (l6.1l): 
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u 


U, 




(17.3) 


Hence the momentum integral becomes 



and thus 

U ( x ) = ^ b pU 0 %( x ) (17.5) 

Hence follows, according to equation (17.1a) 



2 as 

dx 


(17.6) 


V 


On the other hand, one had found before for the smooth pipe, equation (16.13) 
again replacing r by 8 and u^^ by TJ o : 


t q = 0.0225 p.o 7A (^) LA 


(17.7) 


By equating equations (17.6) and (17.7) results: 


i I - °- 022 5 « 


■rar 


This is a differential equation for 5(x). The integration yields: 

4 5A 72 / v ^ A 

5 6 ' T °' 0225 (v-J * 


S(x) = 0.37 ( 0 b * b 




(17.8) 
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- 0-3Tx 


»(*) « ^ 


TJ x 


Re T = — 
z v 


(17-9) 


For the turbulent boundary layer the boundary layer thickness, therefore, 
increases with xV5. The corresponding equation for the l amina r flow 
was, according to equation (9.21a), 6 = 5\/vx/cr 0 . 

By substitution of equation (17*8) into equation (17. 5) one obtains 

a _1/5 

W(x) = 0.036 tpU Q x (HSj.) 

_ g 

or for the drag coefficient c^ = V/~ U 0 xb: 


c f = 0.072 (Be ) 


-1/5 


Comparing this result with test results on plates one finds the numerical 
value 0.072 to be somewhat too low. 



valid for 
5 X 105 < Ee < 10 T 


(17.10) 


corresponds better to the measurements. This law holds true only 
7 

for Re^. < 10 , corresponding to the fact that the Blasius pipe resistance 

law and the 1 /7— power law of the velocity distribution, which form the 
basis of this plate drag law, are not valid for large Reynolds numbers. 

This law is represented in figure 87 together with the laminar— flow law 
according to equation (9.19). The initial laminar flow on the front part 
of the plate can be taken into consideration by a subtraction, according 
to Prandtl (reference 73) 5 


’ -1/5 

Of = 0.07h Re x 


1700 

Re 

x 


5 X10 5 < Re x 
<10 T 


(17.11) 


The plate drag law for very large Reynolds numbers can be obtained 
in essenti ally the s ame way by starting from the universal logarithmic 
law for the velocity distribution equation (16.41) which, according to 
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its derivation ^ is valid up to Reynolds numbers of arbitrary magnitude. 
Here the calculation becomes considerably more complicated. The 
development of the calculation is clarified if one first introduces the 
velocity distribution in a general form. We had introduced for the pipe 

flow the dimensionless variables cp = -H and u = — — . The values at 

v* V 

the edge of the boundary layer are to be denoted by the index 0, thus 


Then: 




u = v # cp 



(17.12) 


(17.13) 


y = 



n 


VCpQ 

fiy - rr— ^ 
u o 


(17-14) 


V 


1 

I'rom. the equation t follows, with W according to equation (17.1) 

° b dx 


2 d , 
P * P < 


u(U 0 - u) dy 


Uo 


dT l A 

v 2 - _d_, 

dx dT] 0 


/»8 


u(U c - u) dy 


and according to equation ( 17 . 13 ): 



Vcp o 

u o 



u(U D - u)dq 


dT »o 


dx 


■s 
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U 


= Vcp 


dT] 


A 'Ho 


o d 


'Po 


2 ° dz dr| c 


— ( i _ JL j dTj 


Ho 


= Vcp, 




dz dt) 


dil 


U 


<Pr 


d n 




dz di} 


uo 




(17.15) 


d \>f qj2\ 

In forming the integral J 1 ( < p — SL )di) one mast note the following 

facts: The differentiation with respect to the upper limit gives zero, 

since for t] = t} 0 , cp = <p 0 . In the differentiation of the integrand cp 

is to he regarded as constant, and cp Q as a function of q Q . Therefore 


l^o 


_<L. 


» - g) an 

\ <Po/ 


Itlo 


l/o 


i»n 0 

dri 

2 

cp 

o o 


Thus follows from equation (17.15): 


i’Ho 


(17.16) 


One puts, for simplification. 


dcp o 

^ ‘ 


A^O 


2 

9 in 


Uo 


(17.17) 


and obtains from equation (17.16): 
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to = rf- F ^o) dT >o 

o 


If one assumes this law to "be valid from the leading edge of the plate 
(x = 0), that is, that the flow is turbulent starting from the front, 
the integration gives: 


with 



(17.18) 


f 1 

iin 0 =° 


(17.19) 


Equation (17.18) can also he written so that the Re-number formed with 
the length of run x appears: 


s®! * T~ ■ 1 


with: 



(17.20) 


(17.19a) 


Equation (17.20) giveB the relation between the dimensionless boundary 
layer thickness t| = v*8 /y and the Re-number U Q X/v » 

The drag remains to be calculated. From 



x 

follows, because W = b J* t q dx 


* 


r 
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3T 


^1 


W = b pU 0 v 


— dri 0 = b p U Q vtC^) 


<P, 


Jt W 


=o 


where 






'(no) 


<p. 


_ dTJ 
2 'o 


^0=° 


The drag coefficient c^ = W/^ U 0 ^b s becomes finally 


c „ = 


2V 


=2 


f V ^i) •(T| 1 ) 


(17.21) 


( 17 . 22 ) 


(17.23) 


Hence c^ also turns out to he a function of Equations (17.20) 

and (17.23) together give a parametric representation of c„ as a function 

r v „8 

of Re x , where the parameter is the boundary layer thickness = — — . 


Numerical Results 

In order to arrive at numerical results, one must introduce a 
special function for ip(q) . For the l/7— power law according to equation 
(l6.ll) that is, with cp = Crp/^, one would obtain the drag law according 
to equation (17. 10) . One uses the universal logarithmic velocity distri- 
bution law, equation (l6.lt). 

<p = 2.5 In q + 5-5 

In order to make the carrying out of the integrations more convenient, 
one writes. 


cp = 2.5 In (1 + 9t0 


Then cp becomes, for tj = 0, cp = 0. The adding of the one chan g es 
cp(T)) a little, only for very small r\ and has only little influence 
on the integrals. If one writes the law at first in the. general form 



38 


MCA TM No. 1218 


<p = a In (1 + Tstj ) 


(17.24) 


the calculation of the integrals equations (17-17)# (17-19)# (17-22), 
with z = 1 + bq yields 

r 

F(r)) = a^ (m 2 z — 2 to z + 2 - 

Z * 

<3>( rj ) = ~ ^z Zn 2 z - 4z In z - 2 ]n z + 6z - 6 j !> (17-25) 


iKtj) = z + 1 - 


2(z ~ 1) 
In z 


With the numerical values 


a = 2.493 b * 8-93 


one obtains for the drag law the followi n g table: 


m 

-6 

Re 10 

B: 

0.500 

0-337 

5.65 

1.00 

0.820 

4.75 

2.00 

I.96 

4.05 

3.00 

3-25 

3-71 

5.00 

6.10 

3-34 

12.0 

17-7 

2.81 

20.0 

32.5 

2-57 

50.0 

96.5 

2.20 

100.0 

217-5 

1.98 

500.0 

1401.0 

1-55 


This table can be replaced by the following interpolation formula: 


°f = 


0.472 


(log Re T ) 


2-58 


ft. 


#* 
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Comparison with test results shows that the agreement Improves if the 
number 0.472 is slightly varied, by putting 



6 9 

10 < Ee I < 10* 


(17.26) 


Prandtl — Schlichting* s Plate Drag Law 

The laminar approach length may again be taken into consideration 
by the same subtraction as before; thus: 


>S 



5 X 10 5 
< Re x < 10 9 


(17.27) 


Whereas the system of formulas equation (17.25) is valid up to Re-numbers 
of arbitrary magnitude, the interpolation formulas, equations (17.26) 

and (17.27) , have the upper limit Re = ICp . However, this limit takes 
care of all Re-numbers occurring in practice. The theoretical formula 
(equation (17.27)) is also plotted in figure 87. Figure 88 gives a 
comparison with test results on plates, wings, and airship bodies. The 
agreement is quite good. 


Very recently this plate drag law has been somewhat improved by 
Schultz-Grunow (reference 89). Until then, the turbulent velocity profile 
measured in the, pipe (l/7-^power law, logarithmic law) had been carried 
over directly to the plate, mainly because accurate velocity distribution 
measurements of the plate boundary layer did not exist. The exact 
measurement of the plate boundary layer showed, however, that the plate 
profile does not completely coincide with the pipe profile. The test 
points show, for large distance from the wall, a slight upward deviation 
from the logarithmic law found for the pipe. Thus the loss of momentum 
on the plate is somewhat smaller than that calculated with the logarithmic 
law. Schultz-Grunow repeated the calculation of the drag law according 
to the formula system given above with the velocity distribution law 
for the plate measured by him. His result is represented by the 
interpolation formula 


°f ~ 


0.427 


(— 0.407 + log Re x ) 


2.64 


10 < Re x < 10" 




( 17 . 28 ) 
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This law is also plotted in figure 87 . The differences from the Prandtl— 
Schlichting law are only slight.* 

The corresponding rotationally— symmetrical problem, that is, the 
turbulent boundary layer on a body of revolution at zero incidence, was 
treated by C. B. Millikan (reference 79)* The l/7— power law of the 
velocity distribution was taken as basis. Application to the general 
case has not yet been made. 


Twelfth Lecture (February 23, 19^2) 
b. The Rough Plate 

The conversion from pipe resistance to the plate drag may be carried 
out for the rough plate in the same manner as described previously for the 
smooth plate. One assumes a plate uniformly covered with the same rough- 
ness k. Since the boundary layer thickness 5 increases from the front 
toward the rear, the ratio k/S which is significant for the drag decreases 
from the. front toward the rear. Behind the initial laminar run, therefore, 
follows at first the region of the fully developed roughness flow; the so- 
called intermediate region follows and farthest toward the rear there is, 
finally, if the plate is long enough, the region of the hydraulically 
smooth flow. These regions are determined by specification of the numerical 
values for the roughness coefficient v#k s /v. In order to obtain the drag 

of the rough plate, one must perform the conversion from pipe flow to plate 
flow for each of these three regions individually. This calculation was 
carried out by Prandtl and Schlichting (reference 76), based on the results 
of Nikuradse (reference 71) for the pipe tests with sand roughness. For 
this conversion one starts from the universal velocity distribution law of 
the rough pipe according to equation (16.28), the quantity B being 
dependent also on the characteristic roughness value v* k Q /v* according 

to figure 8 ^. The calculation takes basically the same course as described 
in detail for the smooth plate in chapter XVTIa. It is, however, rather 
complicated and will not be reproduced here. One obtains as final result 
for the total drag coefficient of the sand— rough plate a diagram (fig. 89 ) 
which represents the drag coefficient as a function of the Reynolds 
number U 0 2/v with the relative roughness l/k 0 as parameter. Just as 
for the pipe, a given relative roughness 2 /k 0 has a drag increasing 
effect not for all Re-numbers, but only above a certain Re-number. This 
diagram is applicable also for roughnesses other than sand roughness, if 
one uses the equivalent sand roughness. In the diagram (fig. 89 ) the 
square drag law is attained. Just as for the pipe, for every relative 


* The tables pertaining to the plate drag formulas are given in 
table 7> chapter XXII. 
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roughness l/k s provided the Re-number is sufficiently large. The 
interpolation formula 


e f = f I.89 + 1.62 log 


jA-2.5 


(17.29) 


applies to this law. 


c. The Admissible Roughness 

The problem of the admissible roughness of a wall in a flow is 
very important in practice since it concerns the effort that might 
reasonably be expended in smoothing a surface for the purpose of drag 
reduction. Admissible roughness signifies that roughness above which a 
drag increase would occur in the given turbulent friction layer (which, 
therefore, still is in effect hydraulically smooth). The admissible 
relative roughness k g /l decreases with increasing Re-number U 0 l/v 

as one can see from figure 89. It is the point where the particular 
curve l/k Q diverges from the curve of the smooth wall. One finds the 
values for the admissible relative roughness according to the following 
table j they can also be combined into the one formula 

U o ke admlss. ^ 10 2 (17-30) 

V 


U l 
0 

V 

10 5 

10 s 

10 7 

10 8 

10 9 

N admiss. 

io -3 

10 * 

10 -5 

10" 6 

10 -7 


From equation (17.3O) one recognizes that the admissible roughness height 
is by no means a function of the plate length. This fact is significant 
for in stance for the admissible roughness of a wing. Equation (17.3O) 
states that for equal velocity the admissible roughness height is the 
aarna for a full scale wing as for a model wing. Let us assume a numerical 
example : 


Wing: chord l = 2m 


Telocity U 0 = 3 00 km/h = 83 m/sec 
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From equation (17. 30) results an admissible magnitude of 
roughness k a = 0.02 mm. This degree of smoothness is not always attained 
hy the wing surfaces manufactured in practice, so that the latter have a 
certain roughness drag. In the considerations just made one deals with 
an increase of the friction drag in an a priori turbulent friction layer. 

However, the roughness may also change the drag by disturbing the 
laminar friction layer to such an extent that the point of laminar/ 
turbulent transition is shifted toward the front. Thereby the drag can 
be increased or reduced according to the shape of the body. The drag is 
increased by this displacement of the transition point if the body in 
question has predominant friction drag (for instance wing profile). The 
drag might be reduced, circumstances permitting, for a body with pre- 
dominant pressure drag (for instance, the circular cylinder). One calls 
the roughness height which causes the transition the "critical roughness 
height”. According to Japanese measurements (reference 77) this critical 
roughness height for the laminar friction layer is given by 


v *^crit 

V 


15 


(17.31) 


A numerical example follows: 

Assume, as prescribed before, a 

wing l = 2m 

U Q = 300 km/h = 83 m/sec 

then Ee = U 0 Z/v = 10^. Consider the point of the wing x = 0.12, thus 
Be x = U Q x/V = 10 . Up to this point the boundary layer might remain 
laminar under the effect of the pressure drop. The wall shearing stress 
for the laminar boundary layer is according to equation (9*17) 

^ = 0.332 U Q 2 = 0.332 £222 x: 2.29 

P K u o x 1C)3 sec 


hence : 

v* = V t o/p = 1*52 m/sec 
and according to equation (17.3I) 


^crit 



-15- I icr^ = 0.14 mm 
1.52 7 
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The critical roughness height causing the transition is, therefore, about 
ten times as high as the roughness height admissible in the turbulent 
friction layer. The laminar friction layer therefore n tolerates M a 
greater roughness than the turbulent one. 

The following can be said about the influence of the roughness on 
the form drag: S h a r p-edged bodies are indifferent to surface roughnesses 

because for them the transition point is fixed by the edges, as for 
instance for the plate normal to the flow. Short cur Ted bodies, on the 
other hand, as for instance the circular cylinder, are sensitiTe. For 
the circular cylinder the critical Reynolds number, for which the known 
large pressure drag reduction occurs, is largely dependent on the 
roughness. With increasing relative roughness k /R (R = radius of the 
circular cylinder) Ee crii - decreases. According to British measurements 

(reference 90 ) the drag curves for a circular cylinder with different 
relative roughnesses have a course as indicated in figure 90. The 
■boundary layer is so disturbed by the roughness that the laminar /turbulent 
transition occurs for a considerably smaller Re-number than for the 
smooth cylinder. The roughness has here the same effect as Prandtl’s 
trip wire, that is, in a certain region of Re-numbers it decreases the 
drag. It is true, however, that the supercritical drag coefficient iB 
then always larger for the rough circular cylinder than for the smooth 
one. 


CHAPTER XVIII. THE TUHBIJLEHT FRICTION LAYER IN 
ACCELERATED AND RETARDED FLOW 


The cases of turbulent friction layer treated so far are relatively 
simple insofar as the velocity outside of the friction layer along the 
wall is constant. Here as for the laminar flow the case of special 
interest is where the velocity of the potential flow is variable along 
the wall (pressure drop and pressure rise). As for laminar flow, the 
form of the boundary layer profile along the wall is variable. In 
practice this case exists for instance for the friction layer on the 
wing, on the turbine blade, and in the diffuser. Of special interest is 
the - question of whether separation of the boundary layer occurs and, 
if so, where the separation point is located. The problem consists 
therefore for a prescribed potential flow in following the turbulent 
friction layer by calculation. The calculation of the turbulent friction 
drag is of importance. The corresponding problem for the laminar friction 
layer was solved by the Pohlhausen method (chapter X). 

For the turbulent friction layer the method of Gruschwitz 
(reference 78) proved best. Gruschwitz makes the asumption that the 
velocity profiles of the turbulent boundary layer for pressure drop and 
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rise can be represented as a one-parameter family, if one plots u/t T 
against y/d« d signifies the momentum thickness which is, according 
to equation (6.32), defined by: 


U 2 ^ = 


u(U - u) dy 


(18.1) 


As form parameter one selects 


1 



(18.2) 


u(d) denoting the velocity in the friction layer at the distance from 
the wall y = d. That rj actually is a serviceable form parameter can 
be recognized from figure 91 where a family of turbulentboundary layer 
profiles is plotted according to Gruschwitz. Gruschwitz found from his 
measurements that the turbulent separation point 1 b given by 


T] = 0.8 (Separation) 


(18.3) 


The form parameter tj is analogous to the Pohlhausen— parameter X of 
the laminar friction layer. However, a considerable difference exists 
between tj and X,: for the laminar friction layer an analytical relation 

exists between X and the pressure gradient and the boundary layer 
thickness, namely according to equation ( 10.41) 


X = 


V 


dH 

dx 


(18.4) 


Such a relation is thus far lacking for the turbulent boundary layer, 
since one does not yet possess an analytical expression for the turbulent 
velocity profiles*. One needs therefore an empirical equivalent for 
equation (18 . 4 ) . 

For the special case of the turbulent' friction layer without 

pressure gradient where the l/T— power law u/CJ =» (y/S ) applies for 
the velocity profile, one finds from equations (18.1) and (l8.2) 


* Compare, however, chapter XXIIb, where under certain assumptions 
such an analytical connection is indicated. 
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| = X , . 0.487 


(18.5) 


Since in the case of the turbulent boundary layer, an analytical 
expression for the Telocity distribution does not exist, the calculation 
is limited to the determination of the -four characteristics of the friction 
layer: form parameter q, wall shearing stress t q , displacement 

thickness 8*, momentum thickness d . Four equations are required for 
their calculation. 

As for the laminar boundary layer, the momentum theorem yields 
the first equation; the momentum theorem may, according to aquation (10. 36) 
be written in the form: 



(18.6) 


The second equation is yielded by the function 



(n) (18.7) 


obtained by Gruschwitz by evaluation of the measured velocity profiles 
(fig. 92), and regarded as generally valid. It can be derived also by 
calculation from the form of the velocity profile (compare appendix 
chapter XXII ) and yields: 


t\ = 1 - 


E - 1 
E(E + 1) 


H— 1 


(18.8) 


The third equation is empirically derived by Gruschwitz from his 
measurements. He considers that the energy variation of a particle moving 
parallel to the wall at the distance v = d is a function of u, d, V. 

Dimension considerations suggest the following relation: 

1 = F(q, Re) (18.9) 

q dx 

q = £ U 2 and S]_ = P + f u 2 ^ signify the total pressure in the layer 
y = d. The evaluation of the test results showed that a dependence on 
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the Re-number is practically non-existent, and that one can represent 
equation (18.9) in the following manner: 

t fSk - 0.00894 T) - 0.00461 (18.10) 

<1 dx 


Furthermore, the identity 


So - g x = P + | U 2 - 


P c 
2 U d 


= flj2 1- 



- <1 T] 


is valid. One puts 


9 n = S 


(18.11) 


and has therefore 


^1 _ dS 

dx dx 

Now equation (18.10) can be written: 


( 18 . 12 ) 


= - 0.00894 g + 0.00461 q 


(III) (18.13) 


The fourth equation is still missing and is replaced by the following 
estimation of t q : According to the calculations for the plate in 

longitudinal flow, equation (17. 7) was: 


T 


o 



0.0225 


us^-A 


-1/4 

= 0.0225 (Be B ) 


(18.14) 


If one takes into ‘consideration that for the l/7— power law of the velocity 
distribution: 
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5* 



■6 = -2- S 
72 


one can write equation (18.14) also: 


— % = 0.01338 

pir 


(R® 8 *) 


-iA 


= 0.01256 (Ee^) 


-lA 


(IT) (18.15) 


For calculation of ■S and q (£ = q tj, respectively) one must now solve 
the following system of equations: 


— + 0.00894 4 = 0.0461 ^ 

+ A + S'N i M = - 1 ° 

<3x \ V <1 dx pU 2 


(18.16) 


q = £ TJ 2 is a given function of x; H and — % are given functions 
2 pU^ 

of tj = £/q or -3, respectively. This system of equations is to he 

solved downstream from the transition point. 

\ 

Initial values: As initial value for i3 one takes the value from 

the laminar friction layer at the transition point: 


^oturh . 


olam. 


(18.17) 


This is based on the consideration that the loss of momentum does not 
vary at the transition point since it gives the drag. The initial value 
of T] is somewhat arbitrary. Gruschwitz takes 

n D = 0.1 

and states that a different choice has little influence on the result. 

With these initial values the system of equations (18.16) may be 
solved graphically, according to a method of Czuber (compare appendix. 
Chapter XXII, where an example Is given) . A first approximation for £ 
is obtained by firBt solving the second equation with constant values 

for r 0 /pu 2 and Hj 
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T n 

= 0.002s H = 1.5 _ (18.18) 

pU 


are appropriate. Thereby the second equation is a differential equation 
of the first order for d. This first approximation is then 

substituted into the first equation, which then becomes a differential 
equation of the first order for £(x)j let its solution be denoted 
by f^(x). Thus one has also a first approximation for tj: T]^(x). 

With T]^(x) one determines the course of H(tj) according to figure 92 
and is now able to improve according to equation (18.15). These 

values of both H and t q are now inserted in the second equation, and 
a second approximation dp(x) is obtained. By substitution of d 2 ( x ) 
into the first equation one obtains the second approximation £ 2 (x), etc. 

The method converges so well that the answer is essentially attained in 
the. second approximation. 

The separation point Is given by 

T] = 0.8 

Incidental to the boundary layer calculation one obtains the following 
characteristic values of the friction layer as functions of the arc 
length x : 


d(x), 5*(x), T}(x), T 0 (x). 

The boundary layer calculation for the profile J 015* c & = 0 Is given 

as example in figure 93 • The transition point was assumed at the velocity 
maximum. The calculation of the laminar boundary layer for the same case 
was indicated in chapter XII. The details of this example are compiled 
in the appendix, chapter XXII. 

It should be mentioned that the calculation for the turbulent 
boundary layer must be performed anew for every Re-number U 0 2/v, whereas 

only one calculation was necessary for the laminar boundary layer. The 
reasons are, first, that the transition point travels with the Re-number, 
and second, that the initial value of d/t varies with Re, since for 

d .fO" 

the laminar boundary layer \ a t the transition point Is fixed. 

It must be noted that the values obtained for t q become incorrect 
in the neighborhood of the separation point: At the separation point t q 

must equal zero, whereas equation (18.15) gives everywhere t q j 0. 
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B oundary Layer Without Pressure Gradient 

Tn this case q(x) = Constant. Equation (18.13) can "be written: 

| = - 0.0089it-Ti + 0.00461 (18.19) 


or, Because q(x) 


constant; 


d£ , d (^ l) = 

dx dr 



(18.20) 


Thus equation (18.19) "becomes: 

d Ji = -. 0.00894n + 0.00461 


(18.21) 


A solution of this equation is: 

1 * - °' 516 (18 ' 22) 


Since at the beginning of the turbulent friction layer T] 1 b s mal ler 
t.bpn this value (transition point T] = 0.1) and since according to 
equation (18.21) dn/dr >0, q must in this case approach the value 
q = O.516 asymptotically from below. For the velocity profile of the 
l/7— power law, q = 0.487 (compare equation (18.15))* The profile 
attained asymptotically for uniform pressure (p = constant) therefore 
almost agrees with the l/7— power law that was previously applied to the 
plate in longitudinal flow. 

A great many boundary layer calculations according to this method 
are performed in the dissertation by Pretsch (reference 80) . 


Thirteenth Lecture (March 2, 1942) 


CHAPTER XIX. FREE TURBULENCE 
a. General Remarks; Estimations 


After considering so far almost exclusively the turbulent flow 
along solid walls, we shall now treat a few cases of the so-called free 
turbulence. By that one understands turbulent flows where no solid walls 
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are present. Examples are the spreading of a jet and its mix ing with the 
surrounding fluid at rest; or the wake flow "behind a "body towed through 
the fluid at rest (fig. 9*0 • Qualitatively these turbulent flows take a 
course similar to that for the laminar case (compare chapter IX) ; quanti- 
tatively, however, considerable differences exist, since the turbulent 
friction is very much larger than the laminar friction. In a certain way, 
the cases of free turbulence are, with respect to calculations, simpler 
than turbulent flows along a wall, since the l amin ar sublayer is not 
present and the laminar friction as compared with the turbulent one can 
therefore be neglected for the entire flow domain. The free turbulence 
may be treated satisfactorily with Prandtl’s concept of the turbulent 
shearing stress according to equation (15*9): 


T 


p r 


Su 

by 


&u 

by 


(19.1) 


the mixing length 2 being assumed a pure position function. The 
turbulent friction has the effect of making the Jet width increase and 
the velocity at its center decrease with increasing distance along the 
Jet. 


We now perform rough calculations, according to Prandtl (reference 2, 
Part I), for a few cases of free turbulence which give information about 
the laws governing the Increase of width and the decrease of "depth" with 
the distance x. 

It has proved useful for such turbulent jet problems to set the 
mixing length 2 proportional to the 'Jet width b; 


1 « 
b 


P = constant 


(19.2) 


Furthermore, the following rule has held true: The increase of the 

width b of* the mixing zone with time is proportional to the fluctuation 
of the transverse velocity v* : 


Db 

Dt 


~ yt 


(19.3) 


D/Dt signifies the substantial derivative; thus: — = 

According to our previous estimation, equation (15.5): 
Therefore : 



v* = 2 
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Bb , du 

— ~ l — 

Bt dy 


Furthermore, the mean value of 


8u 

¥ 


equals approximately: 


^ = number ^!S£ 
dy b 


(19.4) 


and thus: 


Bb 

Bt 


= number 


i V = nuEiber 0 %ax 


(19-5) 


Jet (Plane and Circular) 


We shall estimate, by means of these relations, how the width 
increases with the distance x and the velocity at the center decreases. 
At first, for the circular as well as for the plane Jet: 


Bb 

Bt 


= number u 


db 


max dx 


(19.6) 


It follows, by comparison with equation (19.5): 


— - = number r- = number 
dx b 

b = number x + constant 


If the origin of coordinates is suitably selected (it need not coincide 
with the orifice) one has therefore: 


(plane and circular Jet) (19.7) 


The relation between and x is obtained from the momentum 

theorem. Since the pressure is constant, in the x-direction, the x— 
momentum must be independent of x, thus: 


b = number x 


J = p / d ? = constant 
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whence follows for the circular jet: 


or: 


p p 

J = number p u b 
max 


Vi = nuniber b U “ 


if 


and because of equation (19.7): 



(circular jet) (19*8) 


For the plane Jet, if J* signifies the momentum per unit length of the 
jet: 


J* = number p b 


w = nu2iber — Vt 


1_ 1 /J 7 


and because of equation (19.7): 



(plane Jet) (19*9) 


Wake (plane and circular) 

The calculation for the wake is somewhat different, since the 
momentum which gives directly the drag of the body must be calculated in 
a slightly different way. The momentum integral is now (compare 
equation (9*^0)): 


W = J = 



- u) d F 


(19.10) 


At large distance from the body u* = U 0 — u is small compared with F 0 
(fig, 95) bo that u* « U 0 and 
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u(U 0 - u) = (U Q - u* ) .u* S TT 0 u* 


Thus one obtains for the circular wake: 


§ <V F U 0 2 ~ p XJ Q u* * b2 


c F 


u* w 


IT ,2 

o it b 


Instead of equation (19.6) now applies: 


Db _ n- db 
Dt “ 0 dx 


and instead of equation (19.5)1 


Db , l , 

Dt = nUBiber b u 


Equating of equations (19-12) and (19*13) gives: 


= g U* 

0 dx b 


By comparison with equation (19.II) one obtains: 


b 2 F 

dx n 




F x 


(Wake Circular) 


By insertion in equation (i9.ll) results: 



(19*11) 


(19*12) 


(19.13) 


(19.1k) 


(19.15) 


(Wake circular ) 


(19.16) 
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For the plane wake "behind a long rod., wing, or the like with the 

o 2 

diameter d and the length L, W = c, £■ U_ Ld and 

w 2 u 

¥ = p U Q u» h L 


and hence: 


t 0 d 

Si 


u. 


2 b 


(19.17) 


and from this in combination with equation (19*14): 


2 b 


db 

dx 




~ \|? * 


(wake plane) 


(19-18) 


By substitution in equation (19.17) results: 



(Wake plane) 


(19.19) 


3 

Thus, for the circular wake, the width increases with \/l cl 

- 2/3 V 

the velocity decreases with x ; for the plane wake, the width 


and 


increases with Vx*, and the velocity decreases with x 


rVa. 


The power laws for the width and the velocity at the center are 
compiled once more in the following table. The corresponding laminar 
cases which were partially treated in chapter IX, are included. More- 
over, the case of the free jet boundary is given, that is the mixing of 
a homogeneous air flow with the adjoining air at rest. (Compare 
figure 97 . ) 
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Laminar 

Turbulent 

Case 

Width 

b 

Velocity at Center 
u 

max 
or u’ 

respectively 

Width 

b 

» 

Velocity at Center 
u 

max 
or u’ 

respectively 

Plane Jet 

.00 

CVJ 

H 

x -l/3 

x 

i - 1 / 2 

Circular Jet 

X 

X - ■ ^ 

X 

x-1 

Plane wake 

zl/2 

x-1/2 

X 1 / 2 

X - 1 / 2 

Circular wake 

x 1 / 2 

X -1 

xl/3 

x -2 / 3 

Free Jet boundary 

x 1 / 2 

x° 

X 

x° 


POWER LAWS PCJR TEE INCREASE WITE WIDTH AND THE DECREASE CP VELOCITY" 

whe the distance z 


Eor a few of the oases treated here the Telocity distribution will 
be calculated explicitly below. The calculation on the basis, of the 
Prandtl mixing length theorem was performed for the free Jet boundary, 
the plane Jet, and' the circular Jet by W. Tollmien (reference 8l), for 
the plane wake by H. Schlichting (reference 82 ) and for the circular 
wake by L. M. Swain (reference 83). 


The equations of motion for the plane stationary case are, according 
to equation ( lE . 1^ ) , if the laminar friction terms are completely 
neglected: 


u 




1 dP + 1 Scr z + 1 ^ T xy 
pSc p chc p Sy 


A 


u *I +v *I ss --I& + ± + I 

cbc 8y p $y p 8x p dy 


> (19-20) 


Sv 

3E + 35F 


0 


J 
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b . The Plane Wake Flow 

We shall now calculate the Telocity distribution for the plane wake 
flow. A cylindrical "body of diameter d and span h is considered. 
Further, let 


u* = U — u 
o 


(19.21) 


he the wake Telocity. One applies the momentum theorem to a control 
area according to figure 95, the rear boundary B C of which lies at 
such a large distance from the body that the static pressure there haB 
the undisturbed Talue. As shown in detail in chapter IX, equation (9.40), 
one obtains: 


W = h p / u (TT 0 - u) dy • 
BO 

W = h p J (U 0 — u* ) u* dy 

BC 


(19.22) 


For large distances behind the body u* « U 0 so that one may 

2 

approximately neglect the term u* in comparison with U 0 u* in 
equation (19.22). Hence : 


x oo; W = h p TT 0 / u* dy 

BO 


(19.22a) 


Since, on the other hand W = c.. d h TJ , there becomes: 

w c. O 

l+b 


u* dy = — c d U 

^ 2 w 0 


'-b 


(19.23) 


This problem can be treated with the boundary layer differential equations; 
they read according to equation (19.20) with the Prandtl expression for the 
turbulent shearing stress according to equation (15-9)* with P a x 

neglected: 
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O 

8u + T 3u _ 1 Sr _ 2^2 8u 8 u 
Sx Sy p Sy Sy £y 2 

& + = o 

5x 8y 




> 


J 


For the mir ing length one puts, according to equation (19.2): 


(19-24) 


2 = 3 b 


(19-25) 


Further : 


u* = u^ + Ug + - - . * 


(19-26) 


n 



(19-27) 


For the wake -velocity u-j_ and the width h the power laws for the 
decrease and increase, respectively, with x were already found in 
equations (19.18) and (19.19). One therefore writes: 



(19-28) 

(19-29) 


According to equations (19-27) and (19-29) 


StJ 1 Srj 1 TJ 

^ B( Cw d x) 1 / 2 2 x 


The estimation of the terms in equation (19-24) with respect to their 
order of magnitude in x gives: 


Su —3/2 Sv —3/2 Su —1 S 2 u ^ —3/2 

— ~ x ’ — ~ x ' : — ~ x : ~ x ' : 


8x 


ay 


8y 


X •" V ~ X 


-1 


ay 


*The terms u 2 - - - signify additional terms of higher approximation, 
which disappear according to a higher power of x than does u-j_. 
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-s _2 _o M 

Hence the term v x whereas the largest terms ~ x . Thus 

07 

equation (19.2k) is simplified to: 


du, ? du_ 8 u. 

- U — - = 21 x 


0 Sx 


* * 8 


(19.30) 


The neglected terms are taken into account only in the next approximation. 
The further calculation gives: 


du 

dx 


u ° U <0 

\ v / 


- 1/2 


_ i 1 f 1 i_ f 

2 x 2x 


= tr (■^\ l/2 t i 

^ 0 lo„dJ _ , 1/2 


W / 


B 


(°v d x ) 


A . 


_i = u ^V l/ 2 f” i_ 

Sy 2 \°W d / B 2 C w d X 


2 2 


2 du-^ 8 u^ 


2 / tr 1 ff" 


\ C W 


B 


(V 3 *) 


1 J 2 


After Insertion in equation (19.30) and eliminating the factor 

21/, \~l/2 


0 


Kv 1 ) 


the following differential equation results for f (t) ) 


- (f + n f») = — f»f" 
2 v 1 / B 


(19.31) 


The boundary conditions are: 


y = b: 
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That is: 


il = 1: f = f* = 0 (19.32) 

The differential— equation (19.3O) may immediately he integrated once 
and gives: 

2 

— f ti = — f’^ + Constant 
2 B 


Because of the boundary conditions, the integration constant must equal 
zero; thus: 




This may he integrated in closed form: 


f l/2 ^1/2 = \ 2P_ Af 
B dii 


df , B yl - . 

aw ^ 


2f 


= IS 3/2 + 


20 


2 3 


Constant 




3 V2P 


Because f = 0 for 7] = 


1j C = — i \l~~ and hence: 

3 Y2P 2 


f ^ (1 - ^ 

y 2p 


(19.33) 
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The condition f* = 0 for 73 = 1 is simultaneously satisfied according 
to equation (19*33) • In f", that is , a singularity results at the 

3y 2 ' 


center (ij = 0) and on the edge. For tj = 0, f M = 00 ; the velocity 
profile there has zero radius of curvature. At the edge there exists 
a discontinuity in curvature. In contrast to the laminar boundary 
layer solutions, where the velocity asymptotically approaches the value 
of the potential flow, one obtains here velocity profiles which adjoin 
the potential flow at a finite distance from the center. 


The constant B remains to be determined: 


>+b 


/-b 


Hb I • W | . JL. 

u» dy = *2 u* dy = 2 j u x dy = 2U 0 o^d B j f(rj) drj 

Jo 


0 


From equation (19*33) one finds: 2 ^1 — dTj = 


Iff 


and hence: 


2 / f (ii) dT] = 


B 

20P £ 


and, by comparison with equation (19*23): 


2 



u* 


dy = Oyd 


t>2 

i-T- IT 

2 o 


20P 




or: 

2 

— —r = h B = \flO P (19*34) 

20p 2 2 

Thus the final result for the width of the wake and the velocity 
distribution from equations (19.28) and (19.29) is: 


^Reviewer* s note: Integrating from —1 to +1, as was done in the 

original German version, results in an imaginary term, which was avoided 
in the translation by integrating from 0 to +1 and doubling "the result. 
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a 


b = VlO p (c^d *) 1//2 

^ 1 , m /_ z _\ -1 / 2 r , 

u o " IBP (,o w dJ 



( 19 - 35 ) 


The constant P = 2 /b is the only empirical constant of this theory; it 
must he determined from the measurements. 

Comparison with the tests of Schlichting (reference 82 ) shows tha t 
the two power laws (equations (19.28) and (19.29)) are well satisfied, 
and also bhat the form of the velocity distribution shows good agreement 
with equation ( 19 * 35 )j figure 96 . The constant p is deter min ed as 


P = i = 0.207 

D 


The solution found is a first approximation for large distances; according 
to the measurements it is valid for x/c w d k 50. For smaller distances 

one may calculate additional terms which are proportional to x —1 , x —3 / 2 , 

. . . for the wake velocity in equation (19.26). 

The rotationally— symmetrical wake problem was treated by 
Miss L. M. Swain (reference 83). For the first approximation results 
exactly the same function for the velocity distribution; o nl y the power 
laws for the width b and the velocity at the center u^* are different, 

namely b ~ x 1 / 3 and u^* ~ x -2 / 3 , as already indicated in 
equations (19.15) and (19.16). 


Fourteenth Lecture (March 9, 1942 ) 
c. The Free Jet Boundary 

The plane problem of the mixing of a homogeneous air stream with 
the adjoining air at rest shall also be treated somewhat more accurately 
(fig. 97 ). It is approximately present for instance at the edge of the 
free Jet of a wind tunnel. The problem was solved by Tollmien 
(reference 8l). 

The velocity profiles at various distances x are affine. One 

sets 


u = U 0 f( n ) = U 0 F»(r,) 


( 19 . 36 ) 
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where 


r\ = (h ~ x) 


(19.37) 


and 


^(n) =J f(rj) dq 


Furthermore set 


2 = c x = C b 


(19.38) 


The equation of motion reads: 


dx dy p 3y 


2 j2 du £ g u 

Sy a y a 


( 19 . 39 ) 


One integrates the continuity equation "by the stream function: 

* = f u dy = U Q x f f(T]) di) = U D x F(q) (19.40) 


Then: 


8u 

dx 


U v U >2 U 

^ = h.2 P . 

x ^ 1 ay 2 x 2 


T . - g . _ u o (f - n p') 


(19.40a) 


Substitution into the equation of motion ( 19 . 39 ) gives, after division 
*>y U 0 2 /x; 


FF" + 2c 2 F^F*" = 0 


(19.41) 
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The boundary conditions are: 
at the inner edge: 


t\ = 


at the outer edge: 

il = ti 2 : 


u = TJ : 

F* = 1 

0 

&u „ 
— = 0: 

F" = 0 

oy 

v = 0: 

F = TJ. 

u = 0: 

F* = 0 

^ = 0: 

F” = 0 

3y 


(19. k2) 


Since the boundary points T)-^ and are still free, these five boundary 

conditions can be satisfied by the differential equation of the third order 
equation (19.41). By introduction of the new variable: 


T]* = 




( 19 .^ 3 ) 


the differential equation (19.41) is transformed into (* = differentiation 
with respect to tj*) 


F F" + F" F ,u = 0 


(19.W) 


The solution F" = 0, which gives u = Constant, is eliminated. The 
general solution of the linear differential equation 


F + F*" = 0 


( 19 . 45 ) 


is 


F = e X 1 


* 


Q 

with X signifying the roots of the equation +1=0, thus: 
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A,-] — — 1; 


2,3 


1 i 

2 -+■ I 


V 3 


Hence the general solution is: 


F = C, e 


r*r 



(19-46) 


If, moreover, one measures the q-coordinate from the inner boundary 
point, thus puts: 


T]* = 1}* 



the solution (equation (19-46)) can also be written: 


_TF it 

F = d^ e ^ + d 2 e cos 




al 

2 

+ d^e sin 




From the boundary conditions (equation (19.42)) result for the constants 
the values: 


tj* = O.98I ; tj* = - 2.04 ; » - 3-02 

d x = -0.0062; d 2 = O.987; d 3 = 0.577 


For the width of the mixing region one obtains: 

v3i 


b = x(^T| i - T) 2 ) = x\j2c 2 ^q* 
b = 3-0^[‘ 


-"?) 


/ 2 c 2 x 


The constant c must be determined from experiments. From measurements 
it is found that 


b = 0.255 x 


(19.47) 
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Hence 


and 



0.0845; 


c = 0.0174 


- = 0.0682 

b 


(19-48) 


It ie striking that here the ratio l/b is essentially smaller than for 
the wake. 

The distribution of the velocity components u and v over the 
width of the mixing zone is represented in figure 98 . 

From the second equation of motion one may calculate the pressure 
difference between the air at rest p Q and the homogeneous air stream p^. 

One finds: 


p - p = 0.0048 £ U 2 
i o 2o 


(19.49) 


Thus an excess pressure of one— half percent is present in the jet. For 
the inflow velocity of the entrained air one finds according to 
equation ( 19 • 40a ) : 


= - F(q p ) U 0 = + 0.379 V 20 U < 




and with the measured value of c; 


T_ M = 0.032 TJ q 


(19.49a) 


d. The Plane Jet 

In a similar manner one may also calculate the plane turbulent jet 
flowing from a long narrow slot (compare fig. 94). The laws for the 
increase of the width and the decrease of the center velocity have 
already been given in equations ( 19 « 7 ) and ( 19 . 9 ): b *v x; Uj^ ~ ^=. 

The calculation of the velocity distribution was carried out by Tollmien 
(reference 8l); it leads to a non-linear differential equation of the 
second order the integration of which is rather troublesome. Measurements 
for this case were performed by Forthmann (reference 91). In figure 99 
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the measurements are compared with the theoretical curve. The agre eme nt 
is rather good. Only in the neighborhood of the velocity maximum is 
there a slight systematic deviation. There the theoretical curve is 
more pointed than the measured curve; the ■ theoretical curve, namely, 
again has at the maximum a vanishing radius of curvature. 

According to the Prandtl formula, equation (15.9), the exchange 
■becomes zero at the velocity maximum, whereas actually a small exchange 
is still talcing place. 


e. Connection Between Exchange of Momentum, 

Heat and Material 

In concluding the chapter on turbulent flows I should like to 
make a few remarks about the connection between the turbulent exchange 
and the heat and material transfer in a turbulent flow. 

In the Prandtl theorem equation (15.9) for the apparent turbulent 
stress: 


one can interpret: 



A = p 


du 

du . 


3 r' A 

du 

kg sec 


2 

m 


du 

w 


(19.50) 


as a mixing factor. It has the same dimension as the laminar viscosity p. 
Furthermore, the shearing stress t may be interpreted as a momentum flow: 


_ momentum , „ 

t = = momentum flow 

2 

m sec 


( 19 . 51 ) 


Momentum = mass X velocity = [kg sec] . 

Another effect of the turbulent mixing phenomena, besides the 
increased apparent viscosity by transport of momentum,! s the transport 
of all properties Inherent in flowing matter, as heat, concentration of 
impurities, etc. If this concentration is not uniform, more heat or 
impurity is carried away by the turbulent exchange from the places of 
higher concentration than is brought back from, the places of lower concen- 
tration. Thus there results, on the average, transfer from the places of 
higher bo those of lower concentration. 
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This results, for temperature differences, in a turbulent heat 
transfer; for concentration differences (for instance, of salt), in a 
turbulent diffusion. They can, in analogy to equation (19.5O) be expressed 
as follows; 


Momentum flow = gP°g* = A„ — ^SgntamX 

2 T d.v Vunit massy 

•m na/i ' 0 


m sec 


Heat flow = heat. transport = ^ jL ^ - ‘.heat \ 


m^ sec 


dy \unit mass/ 

Flow of material = ^ter_ial = . ± /material 

2 dy \t 


m sec 


omit mass/ 


The heat content of the unit mass is 

cal m 


c^ = specific heat = 


c p 0 


0 = temperature. 

For chemical or mechanical 


[kg sec degree 

concentrations the concentration of material per unit mass is called the 
concentration c; it is therefore the ratio of two masses and therefore 
dimensionless. Thus the above equations may also be written in the 
following forms; 


T 

Q 


A — 

T ay 



M 



> (19.52) 


J 


The question arises as to whether A T , Aq, are numerically the same 

or different. If the momentum is transported exactly like heat or material 
concentration — Prandtl’s theorem is based on this assumption — it would 
follow that A t = Aq = .A^ and, for instance, the velocity and temperature 

distributions in a turbulent mixing region would have to be equal. 

However, measurements show partially different behavior. 

One has to distinguish between wall turbulence and free turbulence. 
Concerning free turbulence, calculations of G-. I. Taylor (reference 92) 
and measurements of Fage and Falkner (reference 93) showed for the velocity 
and temperature profile of the plane wake flow 



A t 


2 


(free turbulence) 


( 19 . 53 ) 
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The heat exchange is, therefore, larger than the momentum exchange. 
Consequently the temperature profile is wider than the Telocity profile. 
The theory given for that phenomenon by G-. I. Taylor operates with the 
conception that the particles, in their turbulent exchange movements, 
do not maintain their momentum (Prandtl), but their vortex strength 

(PrandtlVs momentum exchange theory ~ Taylor* s vorticity transfer 
theory). However, there are cases not satisfied by the Taylor theory 
(for instance the case of the rotationally symmetrical wake). That the 
heat exchange for free turbulence is considerably larger than the momentum 
exchange is also shown by . experiments of Gran Olsson (reference 88 ) 
concerning the smoothing out of the temperature and velocity distributions 
behind grids of heated rod's. With increasing distance behind the grid 
the temperature differences even out much more rapidly than the differ- 
ences in velocity. 

For wall turbulence the difference between the mixing factors for 
momentum and temperature is smaller. H. Beichardt (reference 87 ) waB 
able to show, from measurements of the temperature distribution in the 
boundary layer on plates in longitudinal flow by Elias (reference 86 ) 
and in pipes by H. Lorenz (reference 14), that here 

A o 

j * = 1.4 to 1.5 (wall turbulence) (19*54) 

T 


Herewith we shall conclude the considerations of free turbulence. 


CHAPTER XX: DETERMINATION OF THE PROFILE DRAG FROM 

THE LOSS OF MOMENTUM 



The method, previously discussed in chapter IX, of determining 
the profile drag from the velocity distribution in the wake is rather 
Important for wind tunnel measurements as well as for flight tests; we 
shall therefore treat it in somewhat more detail. The determination of 
the drag by force measurements is too inaccurate for many cases, in the 
wind tunnel for instance due to the large additional drag of the wire 
suspension; in some cases (flight test) it is altogether impossible. In 
these cases the determination of the drag from the wake offers the only 
serviceable possibility. 

The formula derived before in chapter IX, equation (9.41) for 
determination of the drag from the velocity distribution in the wake Is 
valid only for relatively large distances behind the body. It had been 
assumed that in the rear control plane (test plane) the static pressure 
equals the pressure of the undisturbed flow. However, In practically 
carrying out such tests in the wind tunnel or in flight tests one is 
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forced to approach the "body more closely . Then the static pressure gives 
rise to an additional term in the formula for the drag. For measurements 
close behind the body (for instance.,' for the wing, for x < t) this term 
is of considerable importance, so that it .must be known rather accurately. 

A formula was indicated, first by Betz (reference 84), later by B. M. Jones 
(reference 85 ) which takes this correction into consideration. Although 
at present most measurements are evaluated according to the simpler Jones 
formula, we shall also discuss Betz* formula since its derivation in 
particular is very interesting. 


a. The Method of Betz 

One imagines a control surface surrounding the body as shown in 
figure 100. In the entrance plane I ahead of the body there is flow 
with free— stream total pressure g c , behind the body in plane II, the 
total pressure g 2 <g Q . The lateral boundaries are to lie at so large 

a distance from the body that the flow there is undisturbed. In order 
to satisfy the continuity condition for the control surface the velocity 
ug in plane II must be partially greater than the undisturbed velocity 

U 0 . Consider the plane problemj let the body have the height h. 

Application of the momentum theorem to the control surface gives: 

( 20 . 1 ) 



In order to make this formula useful for test evaluation the integrals 
must be transformed in such a manner that the integrals need to be extended 
only over the "wake”. For the total pressures 


at 


in 


in 


infinity: 

g o 

= P n + | XT ‘ 

O eo 

plane I : 

So 

, P : 

" P 1 + 2 U 1 

plane II: 

g 2 

P 2 

= P 2 + § u : 




> 


( 20 . 2 ) 


/ 


Outside of the wake the total pressure everywhere equals g Q . Hence 
equation (20.1) becomes 



TO 
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W = 



+ 



(20.3) 


Thus the first integral already has the desired form, since the integrand 
differs from zero only within the wake. In order to give the s ame form 
to the second integral, one introduces a hypothetical substitute flow 
Ug' (y) in Plane II which agrees with Ug everywhere outside of the wake, 

hut differs from ug within the wake by the fact that the total pressure 
for Ug’ equals g Q . Thus 


«o = ®2 + I u, 2 2 (20 - lt) 

Since the actual flow u-^, tig satisfies the continuity equation the 

' flow volume across section II for the hypothetical flow u^, u g * is too 

large. It shows a source essentially at the location of the body which 
has the strength 


Q = 



U£» - Ug) dy 


(20.5) 


A source in a frictionless parallel flow experiences a forward thrust 


E = - P U Q Q 


( 20 . 6 ) 


One now again applies the momentum theorem according to equation (20.3) 
for the hypothetical flow with the velocity u^ at the cross section I 

and the velocity u* g at the cross section II . Since g g = g Q and the 
resultant force, according to equation (20.6), equals E, one obtains 


- P TT 0 Q = 




(20.7) 


By subtraction of equation (20.7) from equation (20.3) there results 


W + P U 0 Q = h 


/(« o - * 2 ^ * f/( n *2 2 " "2 8 ) ** 


( 20 . 8 ) 
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or "because of equations (20.5) and (20.6): 




- - 




One may now perform each of these integrations only across the wake, since 
outBide of the wake u*^ = u^. Due to u* 2 ^ — «2^ = (^ u *2 “ U 2)C U *2 + u 2) 
a transformation gives the following formula: 


W = h 

f(e a - g^)ay + f j 

^(u* 2 - ug) 

0**2 + U 2 - OT o) 

dy 


(20.9) 


Betz* Formula 

Tn order to determine ¥ according to this equation, one has to measure 
in the test cross section "behind the body the following values: 

1. Total pressure gg (therewith g Q is the value of gg outside 

of the wake). 

2. Static pressure Pg. 


Furthermore, p = static pressure at infinity. 


Hence one obtains all quantities required for the evaluation of 
equation (20.9)* 


It is useful for the evaluation of wind tunnel tests to introduce 
dimensionless q uant ities. With F = ht as area of reference for the 
drag: 




and hence from equation (20.9): 



( 20 . 10 ) 


For the case in which p^ = P Q = 0 , "k® 8 ^ cro8S section one can 

write this equation, because g Q = q^ 
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This agrees with equation (9.41).* Thus in this case Betz* formula 
changes* as was to he expected* into the previous Biinpie formula. 


b. The Method of Jones 

Later B. M. Jones (reference 85 ) indicated a similar method which 
in its derivation and final formula is somewhat simpler than Betz* method. 

Let cross section II (fig. 101 ) (the test cross section) lie close 
behind the body; there the static pressure pg is still noticeably 

different from the static pressure p Q . Let cross section I be located 
bo far behind the body that the static pressure there equals the undis— 

applies for cross section I according 

(u 0 - u,) dy 1 (20.12) 


back to measurements at cross 
filament is first applied: 


turbed static pressure. Then there 
to equation (9.41)* 



In order to relate the value of u^ 
section II* continuity for a stream 


P dy x = p u 2 dy 2 


(20.13) 


* In chapter IX the total drag of the body (both sides of the 
plates) was designated by 2 W; here the entire drag equals W! 
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Jones makes the further assumption that the flow from cross section II 
to cross section I is without loss* that is, that the total pressure is 
constant along each stream line from H to I: 


g 2 8 1 


(20. lh) 


First, according to equations (20.12) and (20.13): 


W - hp J u 2 (U 0 - u x ) dy. 


(20.15) 


Furthermore : 


So + § n o 2 ‘ B o - «o 


Po + f \ 2 - «! - S 2 


P o 

p + — u <=■ = 
*2 2 2 


«2 


with p =0 
o 


> (20.16) 


and hence 


^2 \ / s 2 ” p 2 


TT 


1o 


Un 


U„ 



(20.17) 


o 2 

From equation (20.12) follows, with ¥ = c^ t h ^ U Q : 


°w 


= 2 


'^2 (, , 12 

V 0 V Vo) d t 


and "because of equation (20.17): 


c w = 2 


|S 2 “ P 2 


Sp \ ^2 

)*-£ 


II 


(20.18) 
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Formula of Jones 


Thus all quantities may be measured in cross section II close to 
the body. This formula is simpler for the evaluation than Betz’ for mul a, 
equation (20.10). 

In the limit, when the static pressure in the test cross section 
becomes p 2 = P Q , this formula, of course, must also transform into the 
simple formula equation (20.11). One obtains for p 2 = P Q = 0 from 
equation ( 20 . 18 ): 



This is in agreement with equation (9.41). 


Fifteenth Lecture (March 16, 1942) 
CHAPTER XXI: ORIGIN OF TOEBTILENCE 

a. General Remarks 


In this section a short summary of the theory of the origin of 
turbulence will be given. The experimental facts concerning laminar/ 
turbulent transition for the pipe flow and for the boundary layer on the 
flat plate have been- discussed in chapter XIII. The position of the tran- 
sition point is extremely important for the drag problem, for instance 
for the friction drag of a wing, since the friction drag depends to a 
great Extent on the position of the transition point. 

The so-called critical Reynolds number determines transition. For 
the pipe (ud/v) cri .j. = 2300 , and for the boundary layer on the plate 

(U 0 x/v) crit = 3 to 5 X 105. However, experimental investigations show 

the value of the critical Reynolds number is very dependent on the 
initial disturbance. The value of ® e cr j_-( : is the higher the smaller the 

Initial disturbance. For the pipe flow the magnitude of the initial 
disturbance is given by the shape of the inlet, for the plate flow by 
the degree of turbulence of the oncoming flow. For the pipe, for instance, 
a critical Reynolds number (ud/v) crit = 40,000 can be attained with very 
special precautionary measures. - 

According to today’s conception regarding the origin of turbulence, 
transition is a stability phenomenon. The laminar flow in itself is a 
solution of th9 Navier-Stokes differential equations up to arbitrarily 
high Reynolds numbers. However, for large Re-numbers the laminar flow 
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"becomes unstable, in the sense that sma1 1 chance disturbances (fluctu- 
ations ,in Telocity) present in the flow increase with time and then alter 
the entire character of the flow. This conception stems from Reynolds 
(reference 101). Accordingly, it ought to "be possible to obtain the 
critical Reynolds number from a stability investigation of the laminar 
flow. 


Theoretical efforts to substantiate these assumptions of Reynolds 
mathematically reach rather far back. Besides Reynolds, Rayleigh 
(reference 102) in particular worked on the problem. These theoretical 
attempts did not meet with success for a long time, that is, no instability 
could be established in the investigated laminar flows. Only very 
recently has success been attained, for certain cases, in the theoretical 
calculation of a critical Reynolds number. 

One ass ume s for the theoretical investigations that upon the basic 
flow which satisfies the’ Navier-Stokes differential equations a disturbance 
motion is superimposed. One then investigates whether the disturbance 
movement vanishes again under the influence of friction or whether it 
increases with time and thus leads to ever growing deviations from the 
basic flow. The following relations will be intorduced for the plane 
case: 


basic flows: 
disturbance movement: 
resultant movement: 


U(3C, y); V(x, y); P(x, y) 
u'(x, y> ; t*(x, y); p'(x, y) 
U + u* . T + v* ; P + p* 


> ( 21 . 1 ) 

J 


P, p* signify pressure. The investigation of the stability of such a 
disturbed movement was carried out essentially according to two different 
methods : 


1. Calculation of the energy of the disturbance movement. 

2. Calculation of the development of the disturbance movement with 

time according to the method of small osc illations. 

I am going to say only very little about the first method since it 
was rather unsuccessful. The second method was considerably more 
successful and will therefore be treated in more detail later. 


The first method waB elaborated mainly by H. A. Lorenz 
(reference 103 )• The following integral expression may be derived for 
the energy balance- of the disturbance movement: 


JL 

Dt 




( 21 . 2 ) 
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In it 1 


= § (u * 2 + v' 2 ) 


signifies the kinetic energy of the disturbance 

movement. The integration is performed over a space which participates 
in the movement of the basic flow and at the boundaries of which the 
velocity equals zero. signifies the substantial derivative. Thus 

one finds on the left side of equation (21.2) the increase with time of 
the energy of the disturbance movement. On the right side. 


M = - 


IT = 


,2 2 57 

u* — + v* — + 

5x 5y 


, ,/So 87 \ 

u'v* — + — ) 
\8y 5x/ 


> (21.3) 





The first integral signifies the energy transfer from the main to the 
secondary movement, the second the dissipation of the energy of the 
secondary movement. If the right side is greater than zero, the intensity 
of the secondary movement increases with time, and the basic flow is thus 
unstable. An assumed disturbance movement u*, v* satisfies merely the 
continuity equation, but no heed is paid to its compatibility with the 
equations of motion. If one could prove that the right side is negative 
for any arbitrary disturbance movement u*, v* this would serve as proof 
of the stability of the basic flow. On the other hand, the instability 
would be proved as soon as the right side is positive for a possible 
disturbance. Unfortunately general investigations in this direction 
are very difficult and have not led to much success. H. A. Lorenz 
(reference 103) treated as an example the Couette— flow (fig. 102), assuming 
an elliptical vortex as a superimposed disturbance movement. He found for 
( 

this case V -S— J, . , = 288, whereas Couette* s measurements for this case 
' v 'krit 
gave the value 1900. 


b. The Method of Small Oscillations 

For the second method (method of small oscillations) the disturbance 
movement is actually calculated, that is, its dependance on the spatial 
coordinates x, y and the time t is developed on the basis of the 
hydrodynamic equations of motion. We shall explain this method of small 
oscillations In the case of a plane flow. In view of the applications 
of this method we shall immediately assume a special basic flow: the 

component IJ, namely, is to be dependent only on y and t and 7 = 0. 

Such basic flows had been previously called ”layer flows”. They exist 
for instance In tunnel flow and pipe flow, approximately, however, also 
In the boundary layer since here ‘the dependence of the velocity component 
U on the longitudinal coordinate x is very much smaller than the 
dependence on the transverse coordinate y. One now assumes a basic flow 
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U(y, t)j TsO; P(x, y) 


(21.4) 




This hasic, flow, hy itself, 
thus 


then satisfies the Navier-Stokes equations. 


su; i sp = v s 2 u 

St + P Sr " Sy 2 



> (21.5) 




A disturbance movement which is also two-dimensional is superimposed 
upon this hasic flow:. 


disturbance motion: u*(z,y,t); Y*(x,y,t); p*(i,y,t) (21.6) 

One then has as the 

resultant motion: u = U + u ’ ; v = 0 + v’; p = P +p* (21.7) 


This resultant motion is required to satisfy the Navier— Stokes differential 
equations and one investigates whether the disturbance motion dies away 
or increases with time. The selection of the initial values of the dis- 
turbance motion is rather arbitrary, but it must of course satisfy the 
c o nt in uity equation. The superimposed disturbances are assumed as ** small", 
in the sense tha t all quadratic terms of the disturbance components are 
neglected relative to the linear terms. According to whether the dis- 
turbance motion fades away or increases with time, the basic flow is called 
stable or unstable. 


By insertion in the Navier-Stokes differential equations (3*18) one 
obtains, neglecting the quadratic terms in the disturbance velocities 

» + + u-Ssi + t« ^ + i SE + i - v + A u 


St St 


Sx 


Sy P Sx P Sx 




St Sx p Sy p Sy 


Su* + Sv* _ q 
S x Sy 


>( 21 . 8 ) 
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If on© now notes the fact that the "basic flow by itself satisfies the 
Navier— Stokes differential equations, equation (21.5), equation (21.8) 
is simplified to: 


8t 


8u’ 

4 _ 

T , M + 

1 

dp* 

8x 

8y 

P 

8x 

8v* + 

u + 

1 

_ 

8t 

Sx 

P 

8y 


8u* 

4 . 

8v* _ 


Sx 

T 

8y 




V A V* 




(21.9) 


The pertinent boundary conditions are: Vanishing of the disturbance 

components u* and v* on the bounding walls. From the system 
equation (21.9) of three equations with three unknown quantities u’ , 
v‘, p‘ one may at first eliminate p‘ by differentiating the first 
equation with respect to y and the second with respect to x and then 
subtracting the second from the first. This gives, with continuity taken 
into consideration: 


> (21.10) 




jSh*- + u &-■ + V* ^ 3 . _ 


^t^y 


= + V 


cbcc^y 

8^u' 


8y 

8-^u’ 


dtox 


u 




,8x 2 8y 


8y 8xdy 


8x c 

8x 3 


In addition to this there is, the continuity equation (21.9). There are 
now two equations with two unknown quantities u*, v*. 


Form of the Disturbance Movement 


For cases where the basic flow predominantly flows in one direction 
as for instance boundary— layer or pipe flow, the disturbance motion is 
assumed to be a wave progressing In the x— direction (= main flow direction), 
the amplitude of which depends solely on y. The continuity equation of 
the disturbance motion may in general be integrated by a disturbance 
function for which the following expression may be used: 
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t (x,y,t) = cp(y)e 


ia(x— ct) 


( 21 . 11 )* 


Where: 


X, = 2jt/a the wave length of the disturbance (a = real) 
c = c r + i 

c r = velocity of wave propagation 

^ = amplification factor; c^ < 0: stable; ^ > 0: unstable 

9(X) - <P r (y) + 1 q> ± (T) - anplltud. of the disturbance »«nt 


From equation (21.11) one obtains for the components of the disturbance 
movement 


u* 


d£ 


9*(y)e 


i<x(x— ct ) 


v* 


chc 


ia(x-ct) 

= - icap(y)e 




/ 


(21.13) 


By substitution into equation (21.10) one obtains the following differ- 
ential equation for the disturbance amplitude cp: 


or 


ia^Dcp" — exp” — cpU” + a cep — Ua cp) = V^cp"" — 2a cp" + a q>) 

(U - c ) (cp" - a^cp) - tf'cp = (cp"” - 2a 2 cp” + Ap) (21.14) 


One introduces dimensionless quantities into this equation by referring 
all velocities to the maximum velocity U of the basic flow (that is 

for the friction layer the potential flow outside of the boundary layer 
and all lengths to a suitable reference length .§) (for instance, for 


*The convenient complex formulation is used here. The real part of 
the flow function, which alone has physical significance, is therefore 






Eefty) = e°^ < 

cp cos 

a ( x " V) 

— cp sin 



L . J 

. 



( 21 . 12 ) 
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the boundary layer flow, the boundary layer thickness) . Furthermore, 
differentiation with respect to the dimensionless quantity y/8 will 
be designated by a prime mark (*). 

Ctne then obtains from equation (21.14) 



(21.15) 


DISTURBANCE DIFFERENTIAL EQUATION 


TJ 5 

where R = -2p. This is the disturbance differential equation for the 

amplitude cp of the disturbance movement. The boundary conditions are, 
for instance, for a boundary layer flow 


y = 0 (wall) : u* 

y = 00 : u* 


v* = 0: cp = cp' = 0 

v* = 0: cp = cp* = o 


\ 


> 


J 


(21.16) 


The stability investigation is an eigenvalue problem of this differential 
equation for the disturbance amplitude qp(y) in the following sense: A 

basic flow U(y) is prescribed which satisfies the Navier-Stokes differ- 
ential equations. Also prescribed is the Eeynolds number E of the basic 
flow and the reciprocal wave length a = 2n/x, of the disturbance movement. 
From the differential equation (21.15) with the boundary conditions 
equation (21.16) the eigenvalue c = c r + i Cj_ is to be determined. The 
sign of the imaginary part of. this characteristic value determines the 
stability of the basic flow. For ^ < 0 the particular flow (U, E) 

is, for the particular disturbance a, stable; for c^ > 0 > unstable. 

The case 0 } = 0 gives the neutrally stable disturbances. One can 
represent the result of the stability calculation for an assumed basic 
flow U(y) in an a, E— plane in such, a manner that a pair of values 
c r , cj: belongs to each point of the a, E— plane. In particular the 

curve c^ = 0 in the a, E— plane separates the stable from the -unstable 
disturbances. It is called the neutral stability curve (fig. 103). In 
view of the test results one expects only stable disturbances to be present 
at small Reynolds numbers for all wave lengths a, unstable disturbances, 
however, for at least a few a at large Eeynolds numbers! ' The tangent 
to the neutral stability curve parallel to the otr-axis gives the critical 
Reynolds number of the respective basic flow (fig. 103). 
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Methods of Solution and. General Properties of the Disturbance 


Differential Equation 

Since the stability limit (e^ = 0) is expected to occur at large 
Reynolds numbers, it suggests itself to suppress the friction terms in 
the general disturbance differential equation and to obtain approximate 
solutions from the so-called frictionless disturbance differential 
equation which reads 


(U - c) (cp" - ct 2 cp) - TJ"<p = 0 


(21.17) 


Only two of the four boundary conditions* equation (21.16), of the complete 
disturbance differential equation can now be satisfied Bince the friction- 
less disturbance differential equation is of the second order. The 
remaining boundary conditions are: 

y = 0: v’ = 0, cp = 0; y = oo; v* = 0: cp = 0 (21.18) 


The cancellation of the friction terms in the disturbance differential 
equation is very serious, because the order of the differential equation 
is thereby lowered from 4 to 2 and thus important properties of the 
general solution of the disturbance differential equation of the fourth 
order possibly are lost. (Compare the previous considerations in 
c ha pter IV concerning the transition from the Navi er— Stokes differential 
equations to potential flow. ) 

An important special solution of equation (21.17) is the one for 
a constant basic flow, U = constant, which is needed for instance for 
the stability investigation of a boundary layer flow as a joining 
solution for an outer potential flow. One obtains from equation (21.17) 
for 


U = constant: 


cp = 



7 


However, due bo the boundary conditions for cp at y = «, the only 
permissible solution is 


cp = e 


-a 7 


(21.19) 


We shall prove at first two general theorems of Rayleigh on the 
neutral and unstable oscillations of the frictionless disturbance 
differential equation. 
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Theorem I: The wave Telocity c r for a Telocity profile with 

U"(y) — 0 mast., for a neutral oscillation (c.^ = 0, c = c r ), equal 
the "basic Telocity at a point so that there exists within the flow a 
point U — c = 0. 

Proof: (indirect) One makes the assumption c > U m (= maximum 

Telocity of the "basic flow). One then forms from equation (21.17) the 
following differential expressions: 

2 U ,t 

L(cp) = 9” - a 9 - 9 = 0 (21.20a) 

U — c 

and 

p_ tt!1 

L(<p) = qp" - a cp cp = 0 (21.20b) 

U - c 


L(cp) signifies the expression obtained from 1,(9), if one inserts 
eTerywhere the conjugate complex quantities. Because of the "boundary 
conditions 


y = 0: 9 = 9 = 0 

y = °°: 9 = qT = o 


One forms further the expression 91,(9) + 9 L(9) and integrates 
it "between the limits y = 0 and y = 00 . The integrals may "be taken up to 
y = oo, since for large y, 9 e -0 ^. Because of equation ( 21 . 20 a, b) 

J-j_ must then be 

noo 


Jn = 




9 L(9) + 9 1(9) dy = 0 




( 21 . 21 ) 


'y=o 


After Insertion of equation (21.20a and b) results, because c” = c, 
p“ 

^9 9" + 9 f, 9 - 2a 2 cp9 — 2 — — — 99^ dy = 0 


U - c J 


y=o 


or 
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* 


J 

1 



+ «p*<p T 



cp’cp’dy 


1 “ 

(■'**) 


qxp dy = 0 


The first term vanishes due to the boundary conditions, hence there 
remains 


0“ r- 


Jl = -2 


qp»q>* 


Uy=o **— 



TJ" 

U - c 


qxp 


dy = 0 


( 21 . 22 ) 


cp'qp’ as well as cpcj) are positive throughout .J if TJ" — 0 and c > TJ^, 
TJ"/U" — c — 0 and hence the integrand in equation (21.22) 1 b positive 
throughout. Thus the integral cannot "become zero. The assumption made 
at the "beginning c > U m therefore leads to a contradiction. 

■< 

For basic flows with TJ — 0, as for instance "boundary— layer flovB 
in a pressure drop, the wave propagation velocity therefore muBt "be smaller 
than TT m for neutral disturbances. Hence a point U — c = 0 exists 
within the flow. This point is a singular point of the frictionless 
disturbance differential equation (21.17) .and plays as such a special role 
for the investigation of this differential equation. The wall distance y 
at which IT — c = 0 is called y = y^. = critical layer. 

This first Rayleigh theorem proved above applies — as shall be noted 
here without proof — in the same manner to flows with U" > 0. 


Sixteenth Lectin’s (March 23, 19^2) 


Theorem II: A necessary condition for the presence of amplified 

oscillations (c^ > 0) is the presence of an inflection point within the 

basic flow (TJ" = 0). 

Proof: (Indirect) According to a ssump tion, ^ ^ Oj thus 

IT _ o ^ 0 for all 'With L(cp) and L(cp) one forms, according to 

equation (21.20a) and (21.20b), a similar expression as before. This 
latter, integrated from y = 0 to y = “, must again give 0, thus 


nco 


J 2 = 



- qp 



dy = 0 


(21.23) 


U 


y=o 
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By substitution according to equation (21.20a and b) results with 
c = c r - i c ± 


or 


J 2 = 


>y=o 


<W>" - W" ~ wf— ^ ~z) 

\JJ - c U - cj 


dy = 0 


= jjpcp* - - 2i 


noo 


U" 


U - c 


W dy ~ 0 (21.24) 


The first term again vanishes because of the boundary conditions. Since cp?p 
is positive throughout and ju — c| ^ 0, the integral can only vanish if 
U < * changes its sign, that is, an inflection point of the velocity profile 
U” = 0 must be present within the flow. It has, therefore, been proved: 

In order to make the presence of amplified oscillations possible, an 
inflection point must exist in the velocity profile of the basic flow, or, 
expressed briefly, such oscillations are possible only for inflection 
point profiles. 

Later on Tollmien (reference 110) proved that the presence of an 
inflection point is not only a necessary but also a sufficient condition 
for the existence of amplified oscillations. Hence the following simple 
statement is valid: Inflection point profiles are unstable. It must be 

mentioned that all these considerations apply in the limiting case E — 
since the proofs were obtained from the frictionless disturbance differ- 
ential equation. 

We know from our previous considerations about the laminar boundary 
layer that Inflection point profiles always exist in the region of pressure 
rise, whereas in the pressure drop region the boundary layer profiles are 
always without an Inflection point, (fig. 104). Hence we recognize that 
the pressure rise or pressure drop 1 b of decisive significance for the 
stability of a boundary layer flow. 

The converse of the theorem just set up is also valid, namely, that 
for E — >oo velocity profiles without inflection point are always stable. 
Bk*om this, however, s one must not conclude that profiles without inflection 
point are stable for all Eeynolds numbers. A closer investigation for 
Beynolds numbers of finite magnitude shows that there profiles without an 
Inflection point also become unstable. One is faced with the peculiar 
fact that the transition from Ee = » to Ee— number of finite magnitude , 
that is, the addition of a small viscosity to a frictionless flow, has a 
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destabilizing effect, whereas one intuitively expects the opposite. As 
later considerations will show in more detail, the typical difference 
between the neutral stability curves of a basic flow with and without 
inflection point appears as represented in figure IO 5 . For the velocity 
profile without an inflection point the lower and the upper branch of the 
neutral curve have, for E — >», the same asymptote a = 0. For the 
velocity profile with inflection point the lower and upper branch of the 
neutral curve have, for E — >», different asymptotes so that for E = 00 
a certain wave length region of unstable disturbances exists. Furthermore 
the critical Eeynolds number is smaller for velocity profiles with an 
inflection point than for those without an inflection point. 

Hence it is to be expected for very large Ee-number, to a first, 
very rough approximation, that the transition point in the boundary layer 
of a body lies at the pressure minimum. Figure 106 shows schematically 
the pressure distribution- for a rather strongly cambered wing profile at 
a sma.l T lift coefficient. The transition point would be expected in this 
case just behind the nose on the pressure side, slightly more toward the 
rear on the suction side. 


Solut i on of the Disturbance Differential Equation 


In order to perform the actual dalculation for the boundary— value 
problem Just formulated, one needs at first a fundamental system 






of the general disturbance differential equation ( 21 . 15 ). 


One imagines the basic flow 
development : 


U(y ) given in the form of a power series 


u(y) = ^o y + -- y 2 


(21.24) 


If one introduces this expression into equation (21.15) and. then wants 
to construct a solution from the complete differential equation which 
satisfies the boundary conditions (equation ( 21 . 16 )), one encounters 
extreme difficulties of calculation, due to the two conditions to be 
.satisfied for y = 00 . In order to obtain any solution at all, one has 
to make various simplifications. The simplifications concern: 

1. The basic flow: Instead of the general Taylor— series 

equation (21.24) one takes only a few terms, thus for instance a linear 
or a quadratic velocity distribution. 

2. The disturbance differential equation: For calculation of ths 

particular solutions the disturbance differential equation is considerably 
simplified. 

Eegarding 1, it should be noted that linear velocity distributions 
frequently have been investigated with respect to stability, as for 
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instance the Couette flow according to figure 102 or a polygonal approxi- 
mation for curved velocity profiles according to figure 107 • This facili- 
tates the calculation due to the fact that then the singular point 
U — c = 0 is avoided in the frictionless disturbance differential 
equation (21.17) for neutral disturbances. However, all investigations 
with linear velocity distributions (references 104, 105, 106) were unsuc- 
cessful with the frictionless as well as with the complete differential 
equation. No critical Reynolds number resulted. When one later took 
for a basis parabolic profiles, these negative results became intelligible. 
One must, therefore, take at least a parabolic distribution as a basis 
for the basic flow. 

Regarding 2, it should be noted that one can provide approximate 
solutions for the solutions of the complete differential equation (21.15) 
from the frictionless differential equation (21.17) Bince the solutions 
are required only for large Re-number R. The frictionlesB differential 
equation however can yield no more than two particular solutions; two 
more have to be calculated, taking the largest friction terms in 
equation (21.15) into consideration. 

The course of the calculation for the particular solutions will be 
briefly indicated. One limits oneself to neutral disturbances, assumes 
a parabolic velocity distribution^ and imagines the latter developed in 
the neighborhood of the critical layer. 


y = y K : U-c-U-c r = 0 

U - c = U* K ^y - + -£- (j ~ (21.25) 

The first pair of solutions cp^, qp^ is then obtained from the friction- 
less disturbance differentia^ equation (21.17) by substitution of 
equation (21.25). According to known theorems about linear differential 
equations with a singular point a linearly independent pair of solutions 
has the form 

U*. (21.26) 

*8 = p s( y - y E ) + 5^ ( y - y z) loe (* - y s) P i( y - y s) 

and Pg are power series with a constant term different from zero. 
The particular solution <p is especially interesting. 
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cp* 



t» for 


j y E 


That is, the u* -component of the disturbance Telocity becomes infi ni tely 
large in the critical layer. This can also be understood directly from 
the frictionless disturbance differential equation (21.17). Accor ding 
to equation ( 21 . 17 ) 


or 


cp” — a^cp 


U" 
U - c 


9 



cp* ~ log fj - y V if TJ” e f 0 


This singular behavior of the solution cp 2 in the critical layer stems 

of course from neglecting the friction. The frictionless differential 
equation here no longer gives a serviceable approximation . In the 
neighborhood of the critical layer the friction must be taken into 
consideration. Moreover, there is another inconvenience connected with 
the cp 2 . For fulfi llme nt of the boundary conditions one requires the 

solution for y - j k > 0 as well as for y - y K < 0. However, for cp 2 

it is at first undeter min ed what branch of the logarithm should be chosen 
at transition from y - y K > 0 to y - y K < 0. This also can be clarified 

only if In the neighborhood of y — y-g-, at least, the large friction terms 

of the complete differential equation ( 21 . 15 ) are taken into consideration. 
The details of the calculation will not be discussed here. The calcu- 
lation leads, as Tollmien (reference 109) has shown, to the result that 
one obtains for the solution cp 2 a so-called transition-substitution in 

the critical layer which appears as follows: 


7 - J K > 0: 

9 2 = P 2 (7 


""k 

(r - *k) 

1 P l(r " 

H 

0 
OQ 

*4“ 

1 

si 

7 ~ 7 E <Q t 

l>> 

V 

ii 

CVI 

e- 

-v s 

UV 
1 + — £ 

( 7 - r x) 

p i( y - p k) 

log | y - y K 





(21.27) 

If one writes, according to this, the complete u'— component, then in the 
neighborhood of y — y-^: 
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7 - J K > °i u 1 



'N 

(as. — pt) 


< 0: u* 



U* 


log 



cos (ax — pt) 


> (21.28) 


+ 


it sin (ax — pt) 





One obtains therefore in the critical layer a phase discontinuity for the 
u*— component. This is retained even in going to the limit, E — > oo . It 
is lost, however, if one neglects the curvature of the "basic flow U" 
or if one operates only with the frictionless differential equation. This 
phase discontinuity is very significant for the development of the motion. 

The loss of the phase discontinuity is the reason that stability investi- 
gations neglecting the curvature U" or operating only with the frictionless 
differential equation remain unsuccessful. 


With this friction correction in the critical layer the pair of 
solutions <p , q>2 * s sufficiently determined. By talcing the friction 

terms in equation (21.15) into consideration, one then obtains a second 
pair of solutions cp^, qp^ which can be represented by Hankel and Bessel 

functions. Of these two solutions cp, tends very strongly towards 


infinity and is therefore not used because of the boundary conditions, 
equation (21.16). <p^ tends, for large y, towards zero. 


The Boundary Value Problem 


The general solution as a linear combination of the four particular 
solutions Is: 


cp = C 1 qp 1 + CgCpg + C 3 9 3 + C^cp 4 


(21.29) 


Let us consider In particular the case where a boundary layer profile is 
investigated with respect to stability. For this case the boundary value 
problem can be somewhat simplified. The previous considerations showed 
that in the disturbance differential equation the friction essentially 
needs to be taken Into consideration only in the neighborhood of the 
critical layer; also, of course, at the wall, because of non-slip. The 
critical layer is always rather close to the wall; hence for y > 5, 
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where U = U m = constant, • one may use the frictionless solution which is 

according to equation (21.19) 9 = e -0 ^. Thus the condition that the 
solution for y = 8 Joins the solution for U = constant is 


‘P'b + a *8 



= 0 


(21.30) 


This mixed boundary condition is therefore to he set up on the outer edge. 
Furthermore, the particular solution 9^ is a priori eliminated in the 

general solution (equation (21. 29)), since it grows, for positive y — y^, 

beyond all limits • 'thus = 0. Hence there re main s for the boundary 

value according to equation (21.16) 


Cm + C 9 + C 9 =0 

1 Y lo 2 Y 2o 3 ^30 


C CD* + C ©’ + C CD* =0 

1 Id 2 2o 3 30 


Vl8 + C 2°26 + C 3*3S = ° 




> 


(21.31) 




A further simplification takes place because of the fact that because of 
the rapid fading away of the solution 9^ on the outer edge y = 8, the 

® already practically equals zero. In the third equation of 
18 


solution 
equation (21. 31) <5 


38 


may therefore be cancelled. The boundary value 


problem actually to be solved is, therefore. 


CD CD CD 

Y lo r 2o y 3o 


ra* 9* m’ 
y lo y 2o y 2o 


$ $ 0 
18 28 U 


= 0 


(21.32) 


This determinant gives the eigenvalue problem indicated above, which 
requires — as has been said before — the solution of the following problem: 
Given 
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1. basic flow TJ(y) 

2. Reynolds number Re = U m &/v 

3. ware length of the disturbance a = 2rt/x 

One seeks from equation (21.32) the pertinent complex eigenvalue 
c = c r + i c ± . Therein gives the velocity of wave propagation and 

c^ the amplification or damping. 

Equation (21.32) may formally be written in the form: 

H(®, c r , c ± , Rj U* o , U" Q , ....) = 0 (21.33) 


where equation (21.33) signifies a complex equation, hence is equivalent 
to two real equations 


If one imagines 
one obtains one 


Hi 

U» , U" , 

0 0 

• 

• 

^ 

11 

0 

E; 

u v u v 

• 

• 

• 

11 

0 

! r 

eliminated 

from these 


(21.34) 


equation between a, R, c^: 



U‘ 0 , u" n . 



(21.35) 


From this equation c-^ can be calculated as a function of a and R. 

The constants U’o, U" 0 , . . . are parameters of the basic flow. Thus, 
if equation (21.35) is assumed solved with respect to c i . 


c i 



U’ 0 , U” D , . 



(21.36) 


Finally one obtains from it, for the neutral disturbances c^ = 0, a 
curve in the a, R— plane, given by the equation 



1 • 


(21.37) 
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This is the sought for neutral stability curve (compare figure 103), which 
separates the unstable from the stable disturbances and. also yields 
the theoretical stability limit, that is, the critical Reynolds n umb er 

Ee crit* 

The performance of the calculation, here only indicated, is analyt- 
ically not possible since the quantities a, c r , R enter into the 

determinant, equation (21.32), in a very complicated manner. One has 
therefore to resort to numerical and graphical methods. The critical 
Reynolds number is very largely dependent on the form of the velocity 
profile of the basic flow, in particular on whether the velocity profile 
of the basic flow has an inflection point, thus on U"(y). 

The critical Reynolds number found from such a calculation gives 
exactly the boundary between stability and instability, hence the first 
occurrence of a neutrally stable disturbance. In comparison with the 
transition point of test results it is therefore to be expected that the 
experimental transition point appears only for larger Reynolds numbers 
where an amplification of the unstable disturbance has already occurred. 


c. Results 

A few results of such stability calculations will be given. The 
completely calculated example concerns the boundary layer on the flat 
plate in longitudinal flow with the laminar velocity profile according 
to Blasius (compare chapter Ha) . In figure 108 the streamline pattern 
of this plate boundary layer with the superimposed disturbance movement 
is given for a special neutral disturbance. Figure 109 shows, for the 
same neutral disturbance, the amplitude distribution and the energy 
balance. Since the disturbance in question is neutral, the energy 
transfer from the main to the secondary movement Is of exactly the same 
magnitude as the dissipation of the energy of the secondary movement. 
Figure 11,0 shows the neutral— stability curve as result of the stability 
calculation according to which the critical Reynolds number Is referred 
to the displacement thickness 6 * of the boundary layer 
(U m 5*/v) crit = 575*' The connection between displacement thickness 5* 

and length of run x is for the laminar boundary layer according to 
equation ( 9 . 21 ) 



Thus a critical Reynolds number formed with the length of run x 


(^m? c /^ ^ crit 
(U m 8*/V) 


crit 


= 1.1 X 10 ^ corresponds to the critical Re-number 
= 575. The critical number for this case observed in tests 
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was 3 to 5 x lo5, . it was explained above that it must he larger t han 
the theoretical number. Furthermore, figure 110 shows t ha t at the 
stability limit the unstable wave lengths are of the order of magn itude 
X = 5$. The unstable disturbances thus have rather long wave lengths. 

This calculation, carried out by Tollmlen (reference 109) for the 
flow without pressure gradient was later applied by Schlichting 
(reference 114, 115) to boundary layer flows with pressure drop and 
pressure rise. The boundary layer profiles with pressure rise and 
pressure drop can be represented in a manner appropriate for the stability 
calculation as a one— parameter family with the form parameter Xp ^ 

according to Pohlhausen's approximate calculation. One then obtains for 
each .profile of this family a neutral-stability curve as indicated in 
figure 111. Hence the critical Reynolds number (U m &*/ V ) cr it ts a 

function of the form parameter Xp ^ according to figure 112. In 
retarded flow (Xp ^ <0) the critical Re-number is smaller than for the 
plate flow (Xp^ = 0), for accelerated flow (X p ^ > 0) it is larger. 

With this result of a universal stability calculation the position of the 
theoretical transition point may be determined conveniently for an 
arbitrary body shape (plane problem) in the following manner: At first, 

one has to calculate for this body the potential flow along the contour, 
furthermore one has to carry out, with this potential flow, a boundary 
layer calculation according to the Pohlhausen method. This calculation 
yields the displacement thickness and the form parameter X p ^ as 

functions of the arc length along the contour, in the form 

“jT = f l^ e ) ^ = f 2^ 


Since in general there exists, accelerated flow at the front of the body 
and retarded flow at the rear Xp^ decreases from the front toward the 

rear. By means of the universal stability calculation according to 
figure 112 one may determine a critical Reynolds number (U m 5*/ V) cr it 

for each point of the contour. The position of the transition point 
for a prescribed Re-number U Q t/v is then given by the condition 



Figure 113 shows, for the example of an elliptic cylinder*, how to find 
the transition point. The curve ( u m 5 */ v ) cr it decreases from the front 

*The boundary layer calculation for this elliptic cylinder was given 
in figure 52 . 



NACA TM No. 1218 


93 


toward the rear; the curve TT m &*/v for a fixed Re-number U Q t/v increases 

from, the front toward the rear. The intersection of the two curves gives 
the theoretical transition point for the respective Re-number U 0 t/v. 

By determining this point of intersection .for various U Q t/v one obtains 
the transition point as a function of U 0 t/v. The result is represented 

in figure 114. The transition point travels with increasing Ee-number 
from the rear toward the front; however, the travel is considerably 
smaller than for the plate in longitudinal flow which is represented in 
figure llh for comparison. Finally figure 115 shows the result of such a 
stability calculation for four different elliptic cylinders in flow 
parallel to the major axis. The shifting of the transition point with 
the Ee-number increases with the slenderness of the cylinder. For the 
circular cylinder the shifting is very slight, which is caused by the 
strongly marked velocity maximum. As a last result, figure 116 shows the 
travel with Ee-number of the transition point on a wing profile for various 
lift coefficients. The profile in question is a symmetrical Joukowsky 
profile with lift coefficients c a = 0 to 1. With increasing angle of 

attack the transition point travels, for fixed Ee-number toward the front 
on the suction side, toward the rear on the pressure side, (compare the 
velocity distirbutions for this profile, given in figure 5h. ) One recog- 
nizes that the shift of the transition point with the lift coefficient 
is essentially determined by the shift of the velocity maximum. 

The last examples have shown that it is possible to calculate 
beforehand the position of the transition point as a function of the 
Ee-number and the lift coefficient for the plane problem of an arbitrary 
body immersed in a flow (particularly a wing) . Regarding the comparison 
with test results it was determined that the experimental transition 
point always lies somewhat further downstream than the theoretical tran- 
sition point. The reason is that between the theoretical and the experi- 
mental transition points lies the region of amplification of the unstable 
disturbances. This amplification also can be calculated on principle 
according to methods similar to those previously described. (Compare 
Schlichting (reference 112) where this was done for the special case of 
the plate in longitudinal flow. ) Presumably one can obtain a still closer 
connection between the theoretical instability point and the experimental 
transition point by applying such calculations to the accelerated and 
retarded flow. 


CHAPTER XXII. CONCERNING- THE CALCULATION OF THE TURBULENT FRICTION 
LAYER ACCORDING TO THE METHOD OF GRUSCHWITZ (REFERENCE 78 ) 
a. Integration of the Differential Equation of the 
Turbulent Boundary Layer 


In order to integrate the system of equations (18.16), one 
first introduces dimensionless variables. One refers the lengths to the 
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wing chord t and the Telocity to the free stream Telocity U 0 , thus: 



d*; 


2 o 





( 22 . 1 ) 


Hence the system of equations (equation (18.16)) may "be written: 


— + 0.00894 — = 0.00461 (- — 

dx* d* \U ' J d* 


d d* 
dx* 





dx* 





> 




( 22 . 2 ) 


First, the second equation is solTed with constant Talues for T 0 /pU^ 
and H, namely 


T 

— ~ = 0.002; H = 1.5 

pu 2 


The first approximation d^*(x*) obtained from that is then substituted 

in the first differential equation. From the latter one obtains a first 
approximation £^*(x*) and from that, according to equation (22.1), 

T]-j-(x*) . With t^Cx*) one determines according to figure 92 the course 

of H(tj) and corrects T q according to equation (18.15). Then me 

obtains from the second differential equation a second approximation 
dg* ( T| ) * etc. 

. For the solution of the differential equations one uses the isocline 
method which can be applied for the present case, according to Czuber, 
in the following manner: Both differential equations haTe the form: 

§ + £(x)y = g(x) (22.3) 

As can be easily shown, this differential equation has the property that 
all line elements on a straight line x = constant radiate from one point. 
The coordinates of this point (= pole) are: 
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£ 


= x + 


1 

f(x) J 


X = fiM 
f(x) 


(22.10 


Thus on© lias only to calculate a sufficient number of these poles and 
can then easily draw the integral curve. 

Figure 93 indicates the result of such a calculation for the profile 
J 015; c a =0. The calculation of the laminar bo undar y layer for the 

same profile was performed in chapter XU, figure 49, table 6. 

Initial values: The transition point was placed somewhat arbitrarily 

at the velocity maximum of the potential flow (x/t = 0.141). It was 
assumed that : 



For the laminar boundary layer was found: 

TJ t 

= I .56 (table 6) 

Hence there results, with 5 */t3 = 2.55; 

^ = 0.141: = 0.611 x 10”3 (table 8) 

The corresponding r] — value was assumed to be 

T] = 0.1 (table 8) 


£.0.141= 22 


Calculation to the second approximation suffices. The result is compiled 
in table 8 and figure 93 • A turbulent separation point does not exist 
since tj remains below 0.8. From the variation of the shearing 
stress t q along the wing chord the drag coefficient of the surface 
friction may be determined: 


t q dx (x = measured along chord) 
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or c v = W/2 t t f U 0 



(22.5) 


The evaluation of tie integral gives 


c v = 0.0090 


b . Connection Between the Form Parameters tj and H = 3 / 8 * 
of the Turbulent Boundary Layer 

According to Pretsch (reference 80) one may alBO represent analyti- 
cally the relation between the form parameters tj = 1 — (u # /u) 1 2 and 
H = 8*/d which was found empirically by Gruschwitz, compare figure 92. 
A power law is set up for the velocity distribution, of the form: 



( 22 . 6 ) 


with n = l/6, 1/7, l/8 . . . , according to the experiments so far. 
Hence results: 


8 * 

T 



Uy/S=0 



(22.7) 




s2n)dz = — L_ 
' n + 1 


1 _ n 

2n + 1 (n + l)(2n + l) 


1 
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From equations (22.7) and. (22.8) follows: 


5 * 

H = z- a 2n + -1 


or 


n 



(22.9) 


From equation (22.8) follows further 


•S H - 1 

ST h(hT +“ l) 


( 22 . 10 ) 


The Gruschwitz form parameter q is defined according to equation (18.2) 
"by: 



With equation (22.6) q becomes: 


n 



( 22 . 11 ) 


Substitution of equation (22.10) into (22.11) gives: 


Tl — 1 

H - 1 

E-l 

q - x 

H(H + 1) 



( 22 . 12 ) 


The connection between H and q calculated according to this equation 
is given in the following table and is also plotted in figure 92. The 
curve calculated according to equation (22.12) almost coincides with the 
curve found empirically by Gruschwitz. 
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H 


1 

0 

1.1 

0.270 

1.2 

0.404 

1.4 

0-573 

1.6 

0.688 

1.8 

0.772 

2.0 

0.833 

2.2 

0.881 

2 .4 

O.916 

2.6 

0.941 

2.8 

0.959 

3-0 

0.972 


Translated hy Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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TABUS vnr. - CALCnLATIOff EXAMPLE TCffi THE TURBULJBfE BOOKDAHY LASBS ACCCKPEK- TO GR13 SC H WITZ (HKFKRB3JC* 78 ) 
EROBTLE J 015; c R =* 0; COMPARE PIC0B3 93- (LAMUAR B0HKDAHT LATER 
PC E THE SAME PBCtFUS; COMPARE PIKDEE 49 AUD TABLE 6.) 


I. ApprcaJjaatlozi 




U. Approximation 




f rW £ 



0.00237 0.611 0.162 0.1 

0.00232 0.70 0.407 0.254 

0.00214 0.87 O.63O 0.406 

0.00211 1.11 0.704 0.478 

0.00200 1.40 0.780 0.521 

0.00191 1.74 0.710 O.553 

O.OOI83 2.17 0.686 0.574 

0.00176 2.65 0.662 0.596 

O.OOI69 3*18 O.635 O.619 

0.00165 3.75 0.611 0.642 

0.00160 4.33 0.593 0.665 

0.00156 4.86 0.577 0.687 

0.00153 5.28 0.565 0.699 

O.OOI52 5.60 0.558 0.708 

0.00151 5.72 0.557 0.715 
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Figure 71.- 


Figure 72. 
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Laminar and turbulent velocity distribution 
in pipe. 
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Fluctuation with time of the velocity of 
turbulent flow at a fixed position. 
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Figure 75.- Laminar and turbulent boundary layer on a flat plate in 

longitudinal flow. 
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Figure 76.- Transfer of momentum by the 
turbulent fluctuation velocity. 
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Figure 'll.- Explanation of the mixing length. 
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Figure 81.- Universal resistance law of the smooth pipe. 
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Figure 82.- Resistance law of the rough pipe. 
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Figure 83.- Resistance law of the rough pipe for fully developed 

roughness flow. 



Figure 84.- The roughness function B as a function of v kg /v. 
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(a) Roughness tunnel. 



(b) Velocity distribution in the section A A. 

Figure 85.- Measurement of the drag of an arbitrary roughness. 
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Figure 86.- Calculation of the turbulent plate drag. 
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Figure 88.- Ibe drag law of the smooth plate; comparison with measurements. 
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figure 89.- The drag law of tile rough plate. 
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Figure 90.- Plane wake flow; comparison of measurement and calculation. 
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Figure 98.- Free jet boundary; distribution of the longitudinal and 

transverse velocity. 
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Figure 100.- Determination of the profile drag according to Betz. 



Figure 102.- The stability investigation of the Couette flow. 



Figure 103.- The neutral stability curve as result of a stability 

investigation (schematic). 



Figure 104.- Basic’ flow without and with inflection point. 
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Figure 105.- Neutral stability curves for velocity profiles without and 

with inflection point (schematic). 



D = Pressure side 
5 = Suction side 
U = Transition point 

Figure 106.- Wing profile with pressure distribution. 
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Figure 107.- Approximation of a velocity- 
profile by a polygon. 



Figure 108.- Streamline pattern and velocity distribution of the plate 
boundary layer for neutrally stable disturbance U(y) = basic flow; 
u(y) = U(y) + u'(x,y,t) = disturbed velocity distribution; 

A = 27r/a = wave length of the disturbance. 
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Figure 109.- Neutral disturbance for the friction layer on the flat plate 
S Amplitude of the disturbance velocity: 

u'(x,y,t) = u x (y) cos (x - P r t) - u 2 (y) sin (x - P r t)- 
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Figure 110.- Neutral stability curve for the friction layer on the flat 

plate in longitudinal flow. 
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Figure ill.- Neutral 
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7 curves for the boundary-layer profiles in 
rated and retarded flow. 
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Figure 113.- Stability calculation for the elliptic cylinder of axis 

ratio 3-i/b^ = 4. 



Figure 114.- Result of the stability calculation for the elliptic cylinder 

of axis ratio a iA>i = 4. 
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A = laminar separation point (independent of R e = 
M = maximum velocity of the potential flow. 


Figure 115.- The position of the instability point as a function of the 
Reynolds number for the elliptic cylinders of axis ratios 
a l/b^ — 1, 2, 4, 8. 
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Symmetrical Joukowsky profile: 1 = 0.15 



A = laminar separation point (independent of R e = -^- ) 
M = maximum velocity of the potential flow 
S = stagnation point 


Figure 116.- The position of the instability point as a function of the 
Reynolds number for a Joukowsky profile for lift coefficients of 
c a = 0 to c a = !. 



